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Symbol 

A 

A 

a = 

a = 

b = 

Ck 

CD = 

Cd 

CF = 

c = 

CP = 

D = 

d = 

E = 

e = 

e = 

F = 

F 

F = 


NOMENCLATURE 

Description 

Component of roller raceway interference at 
aximuth j ; m ( in . ) 

2 

Area, m2 (in2 or ft ) 

Semi-major axis of a contact ellipse; m (in.) 
Bearing fatigue life adjustment factors 
Semi-minor axis of a contact ellipse; m (in.) 
Viscosity constants = 1, 2 and 3 
Crown drop; m (in.) 

Drag coefficient 
Centrifugal force 
Constant; cj = 1, c 2 = -1 
Specific heat - BTU/lb.°F 

Rolling element or general diameter; m(in) 
Cage rail and ring land diameters; m(in ) 
Young's modulus.; psi 
Pitch diameter; in. 

Surface eraissivity 
Friction Force; N(lb) 

Applied force or moment; N(lb) or N-m(in.-lb) 
Shaft and housing fits; m tight (in. tight) 
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Symbol 

G 

g 

o 

H 

H 


h 

h 

J 

J 

K 

K 

K 

K 

L 

L cm 


>-10 


Description 
Gyroscopic moment; N-m(in.-lb) 

Gravitational constant; m/sec^ (in./sec2) 

Heat generation rate; Watts (Btu/hr.) 

Height of the roller end flange contact above the 
inner raceway measured along a line defining the flange 
contact angle; m(in) 

Lubricant film thickness; m(in) 

Convective heat transfer coefficient; Watts/m2°K, (BTU/ 
Hr-ft^°F) 

Mass moment of inertia; m-Kg-sec^, (in-lb-sec^) 

Mechanical to heat energy conversion factor 
Thermal conductivity; watts/m-°K (BTU/hr-f t- °F) 

Ratio, Roller sphere end radius to rearing apex length 
Lubricant mass flow rate proportionality factor 
Integration constant 

Roller total length or a general length; m(in or ft) 
Length from roller small end to the center of mass; 
m ( i n ) 

Bearing fatugue life, at which 90 percent of the bearings 
survive (Hr) or revolutions 
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NOMENCLATURE (CQNTINOEU) 


M 

o 

M 


Symbol Description 

^ = Roller raceway effective contact length; m(in,) 

= Moment; N-m(in-lb) 

= Fluid mass flow rate; Kg/Hr, (Ib/Hr) 

N = Total number of ro 1 1 er/raceway lamina jets, under 

race flow paths 

N = Rotational speed; RPM 

N = Rotational speed factor 1/m^ (l/in^) 

NRe = Reynolds Number 

P = Decimal fraction volume of lubricant/volume of 

bearing cavity 

P = Shaft/inner ring, housing/outer ring, external 

2 

pressures N/m (psi) 

2 

P* = Change in fit pressure due to ring inertia; N/m (psi) 

Pr = Prandtl Number 

Pd = Diametral clearance; m, (in) 

0 = Contact normal load; n (lb) 

Qc = Raceway dynamic capacity; N (lb) 

i^e = Raceway equivalent load N (lb) 

Q = Heat transfer rate; watts (BTU/Hr) 

0* = Heat generated at a given temperature node; Watts 

(BTU/Hr) 
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Symbol 

q = 

R = 

R = 

Rsph = 
r = 

rf 

S = 

s = 

T = 

T 

t = 

U 

u = 

0 = 

u = 

V = 

V = 


NOMENCLATURE (CONTINUED) 

Description 

Load per unit length; N/m Ib/in 
Radius; ra (in) 

Roller/raceway equivalent Radius; ra (in) 

Roller end sphere radius m(in) 

Radius ; m (in) 

Radius from roller axis of rotation to the center of the 
roller end flange contact; m (in) 

Dimensionless speed factor 
y/a 0 s < 1 

Temperature; °K (°F) 

Bearing shaft andhousing mean temperatures; °K ( °F ) 
z/b 0 t < 1 

Radial portion of contact deformation in a roller/ 
raceway contact; ra (in) 

Uniform radial displacement; m (in) 

Dimensionless entrainment velocity 

Mean value of contact surface velocities; m/sec (in/sec) 
Sliding velocity; m/sec, (in/sec) 

Fluid velocity; m/sec (in/sec) 
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Symbol 

W 

X 

Y = 

Z = 

o< = 

^ = 

= 

/ 

■-C 

= 

A 

A< = 

a = 

<3 
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NOMENCLATURE (CONTINUED) 

Description 

Width or thickness; in 

Dimensionless factor to ad just bearing shaftor housing 
shaft rigidity accounting for unequal bearing ring 
and shaft and housing effective widths 

Dimensionless factor to adjust bearing shaft or housing 
seat rigidity accounting for shaft and housing shoulders 
Number of rollers 

Contact angle or roller included angle; degrees or 
radians 

Coefficient of thermal expansion; (44rrp—) 

+- _ Tl degrees or radians 

D Co s°<m 
e 

D 

e 

Lubricant specific weight at 60°F 
Designates change 

Inner ring displacements relative to outer ring; m(in) 

VC. - 1- -a. 

Roller raceway deflection; m(in) 

Difference between ring radial displacements at a 
common boundary due to rotational speed; m(in) 

Exponent for calculating raceway equivalent load 
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Symbol 

€ 

^ = 

? = 
r = 

J\. = 

X = 

X = 
} = 

/ 

e 

If = 

^ = 
^ = 


NOMENCLATURE (CONTINUED) 

Description 

Lubricant scoop or Jet efficiency 
Roller pitch angle degrees or radians 
Absolute viscosity; N-sec/m^ (ib^ - sec/in^) 
Absolute Viscosity; Kg/Hr-m Ib/Hr-ft 
Time; seconds 

Ratio; lubricant film thickness to composit surface 
roughne ss 

Pressure coefficient of viscosity; m^/N (in^/lb^) 

Heat transfer coeii icietit ; watts/'"K . IbTU/Hr i . 
Constants = o ">2 “ ^ 

Poissons Rat i o 
Density; Kg/ra^ (Ib/in^) 

Pressure; N/m^ (psi) 

Surface roughness; ra(in) 

Composite roller/raceway surface roughness 
Curvature sura; l/ra (l/in) 

Shear stress; N/m^ (psi) 

El 1 ipcitivity Factor 

Roller axirauth angle; degrees or radians 
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NOMENCLATURE (CONTINUED) 


Symbol 


Description 


Cp 

JL 

Ui 

0)c 


Ratio of jet travel to roller length 

Ratio of lubricant flow rates through under-race 
flow paths. 

Inner or outer ring rotational speed radians/sec 
Rotational speed relative to cage speed; radians/sec 
Cage rotational speed; radians/sec 


= Roller orbital speed; radians/sec 


Subscripts 

B = Refers to lubricant distribution 


cold = 

CR = 

c = 

d 

F 

(F.Q) = 


f = 

H = 
HI = 


refers to shaft and housing fits at 293®K (68°F) 

Cage Rail 
Refers to cage 

Refers to lubricant hydrodynamic drag 
Refers to lubricant distribution 

Refers to ring radial displacement due to fits and 
uniform rolling element loads. 

Refers to flange 
Refers to housing 

Refers to the housing outer ring equivalent section at 
the roller flange end side of the bearings 


Ratio of lubricant flow rates through under race symbol 
flow path? 
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S ymbo 1 
H2 

Hot 

I 

i 

J 

j 

k 

LL 


LT 


m 


n 


NOMENCLATURE (CONTINUED) 

Description 

Refers to the housing outer ring equivalent section 
at the roller small end side of the bearing 
Refers to shaft and housing fits at operating 
temperature 
Refer to inner ring 

Refers to inner ring, internal pressure or to a temper- 
ature node 
Refers to jet 

refers to a roller azimuth location, or to a temperature 
node 

Refers to a laminum location or to a temperature node 
Refers to ring radial expansion which occurs, starting 
with loose fit, up to the speed at which the fit 
becomes tight 

Refers to ring radial expansion which occurs, starting 
with a loose fit, at the total ring speed 
m 1 refers to the outer raceway 
m 2 refers to the inner raceway 
m 3 Refers to the flange 

n-1 refers to the roller large end side of the bearing 
nz2 refers to the small end side 
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Symbol 

0 

0 

R 

S 

S 

51 

52 

SP 

XL 


T 

T 


TOT 


NOMENCLATURE (CONTINUED) 

Description 

Refers to outer ring or external pressure 

Refers to outer ring 

Refers to roller 

Refers to shaft 
Refers to scoop 

Refers to inner ring shaft equivalent section at 
the roller large end side of the bearing 
Refers to the inner ring shaft equivalent section 
at the roller small end side of the bearing 
Refers to speed related lubricant distribution 
Refers to ring radial expansion which occurs; 
starting with a tight fit, up to the speed at 
which the fit becomes loose 
Refers to total specified ring space 

Refers to ring radial displacement which results from 
changes in temperature from 293°K (68°F) 

Refers to total ring radial displacement which results 
from fit, temperature, rolling element load and 
rotational effects 
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Symbol 


NOMENCLATURE (CONTINUED) 

Description 


X 

y 

z 

n. 

1 

2 

3 

4 

5 


Refers to a distance measured in the x direction 

Refers to a distance measured in the y direction 

Refers to a distance measured in the z direction 

Refers to rotational speed 

Refers to outer raceway, x direction, or roller large 
end side of the bearing 

Refers to inner raceway, y direction, or the roller 

small end side of the bearing 

Refers to flange or z direction 

Refers to the XZ plane 

Refers to XY plane 


XX 
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1 . SUMMARY 


The purpose of this effort was to develop a tool for the 
analysis of high speed tapered roller bearings. The tool 
developed i s S CS [F Computer Program AE71Y001 which has the 
following capabilities. 

1) Bearing loads may be applied and bearing deflections 
determined in three distinct linear directions and two 
angular directions. 

2) Hydrodynamic cage and roller drag friction and elasto- 
hydrodynamic concentrated contact (EHD) friction forces 
are considered. 

3) Cage web/roller loading is determined. 

4) Bearing thermal performance is determined on the basis 
of heat generation, dissipation rates and temperature 
mapping. 

5) Bearing geometry changes which occur from cold to 
operating conditions are determined as affected by thermal 
gradients, shaft and housing fits, ring centrifugal growth, 
and rolling element - raceway normal loads. 

As previously available to S G§ JF, the computer program had 
capabilities 1, 2 and 3. Under Contract NAS3-15700 capabilities 
4 and 5 were added. 

The complete computer program was used to optimize the design 
of a tapered roller bearing having the external dimensions and 
operating conditions listed below. 

Bearing Boundary Dimensions: 

Outer ring O.D. - 2061375 mm (8.125 in) 

Inner ring I.D. - 120.650 mm (4.750 in) 

Total Width - 47.625 mm (1.875 in) 

Operating Conditions: 

Axial Load - 8896-53389 N (200-12000 lb) 

Radial Load - 8896-26678 N (2000-6000 lb) 

Shaft Speed - 12500 RPM 
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Lubricant Data: 

Type-MIL-L-23699 

Input Temperature - 364 °K (195°F) 

Maximum Lubricant Flow Rate - .017 £1.1 — (4.5 gal/min) 

min 3 ' 

Based on the maximum load conditions, the optimum internal 
geometry was developed as follows: 


Outer Raceway Half Angle 
Inner Raceway Half Angle 
Roller Included Angle 
Roller Large End Diameter 
Roller Length 
Roller End Sphere Radius 
Number of Rollers 


23° 57’ 36" 

18° 20’ 

5° 40’ 

23.609mm (0.9295 in) 
27.742mm (1.0922) 
196,850mm (7.75 in) 
20 


The bearing half angles are biased by 2’ 24" such that roller large 
end pinching is promoted. This bias is relieved at operating 
conditions. 


The bearing lubrication scheme was also studies and optimized. 
Recommended lubricant flow rates are 227 Kg/Hr (500 Ib/Hr or about 
1 gal/min) introduced to the inner raceway undercut at the roller 
large end and 227 Kg/Hr via jets, to the inner race impinging at 
the roller small end side of the bearing. The lubricant is to be 
scavenged at both sides of the bearings. 

Under the maximum load conditions, and with the recommended 
lubrication scheme, the predicted lubricant temperatures within 
the bearing are less than 422°K (300°F). 

For the optimized bearing at the maximum load conditions, 
the bearing Lj^q life is 765 hours. This includes a 2.55 
mi 1 1 ip 1 i cat ive factor which reflects practically full EHD films 
in the roller raceway contacts. A material life multiplier of 
1.0 was assumed. Cage speed is predicted to deviate from eppi- 
cylic by less than 6 percent. Under the same conditions the roller 
end flange contact performance parameters are as follows: 

Maximum Hertz stress -3.01 x 10*^ N/m^ (43,750 psi) 

Maximum Sliding speed 22.85 m/sec (900 in/sec) 

Maximum Rolling Speed 49.50 m/sec (1950 in/sec) 

Minimum EHD film thickness 0.965 m x 10“^, (38 mi c ro i n ch e s ) . 
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A "Z" configuration cage which is guided by the outer 
ring I.D. at the roller small end side of the bearing and at the 
inner ring flange at the roller large end side of the bearing 
was selected as optimum because the lubricant flow to and from 
the bearing is less restricted than it would be with either the 
strictly inner or outer ring land riding cage. 
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2. INTRODUCTION 

The purpose of this effort was two fold: 

1. To develop a mathematical model in the form of a computer 
program for the determination of the kinetic and thermal 
performance characteristics of lubricated high speed 
tapered roller bearings. 

2. To use the mathematical model to optimize the internal 
geometry of a tapered roller bearing whose boundary 
dimensions and operating conditions are presented in 
Enclosure 1. 

The basis for the model was embodied in the S C§ IF developed 
Computer program AE71Y001, "The Analysis of Dynamic Performance 
Characteristics of Tapered Cylindrical and Spherical Roller 
Bearings Under Radial Axial and Moment and/or Misalignment 
Loading". Under the sponsorship of NASA, a heat dissipation, 
temperature mapping analysis was developed, and incorporated into 
this roller bearing computer program. In addition, an S CS [F 
developed model which predicts the change in bearing geometry 
from the "off the shelf" state, to operating conditions was 
incorporated into the program. As finally constructed, computer 
program AE71Y001 can analyze total bearing performance when the 
bearing is subjected to general loading, high speed rotation and 
various combinations of jet and through-race lubrication. 

The methods of analysis used are extensions of previous 
work. The bearing load def 1 e ct i on an a ly s i s wa s first developed 
by Jones ref. (1). Harris expanded the load deflection analysis 
to include the effects of hydrodynamic (HD) and el astohydrodynami c 
(EHD) lubrication on rolling element bearing performance, refs. 

(2 and 3). Crecelius and Harris utilized a temperature mapping 
heat dissipation model in conjunction with a bearing dynamics 
model, to optimize the design of an ultra high speed ball bearing, 
ref. (4). The analysis presented herein offers further refinements 
to the previous work by providing better definition of lubricant 
distribution within the bearing cavity, and a more complete 
analysis of the changes in bearing geometry which occur as a bearing 
reaches stea (fystate operating conditions. 

It is expected that these models will be further refined 
as analytical methods improve and comparisons between analytical 
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predictions and experimental results become available. 

Interest in high speed tapered roller bearings has been 
provided by the need for greater load carrying capacity in 
moderately high speed gearbox applications, typically serviced by 
ball bearings. The immediate need for high speed, high torque 
gearboxes has been generated by the helicopter industry. Addi- 
tionally, tapered roller bearings are being evaluated for use 
as raainshaft bearings in small gas turbines. 

Results from several high speed tapered roller bearing 
experimental programs (ref. 5 and 6) have been used in the 
bearing design optimization portion of this effort. It 
became clear as a result of this optimization study, that 
additional experimental and theoretical work would be required to 
develop a complete model defining acceptable roller-end-flange 
contact. This analysis provides various operating characteristics 
realized at that contact such as: 

• Load 

• Hertz Stress 

• Sliding Speed 

• Rolling Speed 

• EHD Film Thickness 

• Lubricant and Surface Ambient Temperatures 

However, a unifying theory describing a successful REFC, in 
terras of the above variables does not exist. 

The bearing performance data includes load/deflection data 
in the form of spring constants for each of the bearing’s five 
degrees of freedom. This data is useful for rotor dynamics 
studies, and its presentation was an explicit requirement of this 
effort. 

Throughout this analysis English Units are used. Where data 
is presented in International Standard Units, it was obtained 
by conversion utilizing ref. (7). 
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3.0 ANALYSIS 


3.1 Mathematical Models 


The mathematical model has four major segments which are 
interdependent but sufficiently distinct to allow them to be 
discussed separately. The four analytical segments and their 
major subsegments are as follows: 

1. Bearing dynamic performance 

a. Applied loading, bearing deflection and inner ring 
force and moment equilibrium. 

b. Lubricant related friction forces and roller and 
cage equilibrium. 

2. Thermal performance 

a. Temperature mapping 

b. Heat dissipation 

c. Lubricant distribution 

3. "Off the shelf" to operating geometry analysis. 

4. Bearing raceway fatigue life 

The size and diversity of this analytical presentation has 
required some designations to be used with several meanings. 
Therefore, in the required instances, their meanings are pre- 
sented immediately subsequent to their use. 

3.1.1 Bearing Dynamic Model 

A major portion of the following presentation is taken 
directly from a publication. Rolling Element Bearing Dynamics 
by Harris and Mindel, ref. (8). That work is not included in its 
entirety as an appendix because it deals with ball, as well as 
roller bearings. In addition. Rolling Element Bearing Dynamics 
presents a more detailed analysis of the cage-roller interaction 
than is used in computer program AE71Y001. In particular, the 
AE71Y001 analysis assumes zero friction at the roller/cage web 
contact . 
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3. 1.1.1 Roller Analysis 


The nature of the tapered roller bearing dynamics problem 
is such that it is most easily understood if the action of a 
rolling element is examined, first, then, the interaction of each 
rolling element with the cage, rings and lubricant will be 
developed. 

3. 1.1. 1.1 Normal Forces and Deformations 


The tapered roller bearing and the individual roller coordinate 
systems are given in Figure 1, 

A roller may be characterized by Figure 2; it rotates about 
an axis inclined to the x-axis and its lateral surface is a 
portion of a curved surface of revolution. This surface is 
arbitrarily shaped, and for the purpose of analysis a roller 
is divided into a number of laminae, each having identical width 
W„j, It is assumed that there are no abrupt pressure changes along 
the width of the contact and thus no shear effects are assumed 
to occur between laminae and each laminum deflects independently 
of its neighbors. 

In the following equations, the sub-scripts take on the 
following meaning: 

j - refers to a particular roller j = 1..Z 

where Z is the total number of rollers in the bearing 

m - refers to a particular raceway 
m = 1 outer raceway 
ra = 2 inner raceway 

k - refers to a particular laminum location k=l...N where N 
is at most 20. 


Lundberg ref. (9) developed the following load-deflection 
relationship for a roller in contact with a raceway such that a 
uniform load per unit length occurs along the length of contact: 


C . = 2 0 m 1 ( ' ? ) ) „ 


• mj 


ITE-i 


r?J- 


3 . 1 . 1 . 1 


Qoaj {l“ ^ 1 4- C m Yno)J 

In terms of load per unit length acting on a particular laminum,/ 
Eq. 3. 1.1.1 becomes: 
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3. 1. 1.2 


whereXmj is the actual length of contact. The roller load is: 


\c--nJ 







3. 1. 1. 3 


f is Poisson's ratio, E is the modulus of elasticity. Both 
roller and raceway materials are assumed to be steel. 


Ym 




3, 1. 1.4 


where D is the roller diameter, e is the pitch d iamet er , s the 
raceway contact angle and c., =+l ^ Cl = -1 . 

3. 1.1. 1.1.1 Roller Alignment 

In a tapered roller bearing, rollers are constrained by the 
raceways and guide flange to move in. a prescribed path and the 
inner and outer raceway contact angles are substantially constant. 
However, a roller axis can pitch through a small angled ; relative 
to the contact angle to achieve a balance of moments about the 
roller centroid in a plane passing simultaneously through the roller 
and bearing axes of rotation. The maximum normal displacement 
owing to roller pitch is equal to i ra j j . The skew angle is 
assumed zero. 

3. 1.1. 1.1.2 Crowning 

When a tapered roller is pressed into a reacway of different 
axial length, stress concentrations develop at the contact 
extremities as shown by Figure 3. To relieve these stresses, 
crowning,' shown by Figure 4, is used. The crown radius is usually 
large compared to the roller diameter, so that the crown drop at 
the roller ends is slight, but sufficient to minimize edge 
loading . 

3. 1.1. 1.1. 3 Contact Deformation with Roller Pitch and Crowning 

At a given azimuth angle^^j (see Figure 1), the contact 
deformation varies as a function of laminum position. Assuming 


8 

RESEARCH LABORATORY SKP INDUSTRIES, INC. 



F-4I8SA RIM 


AL73P010 


a contact is divided into N laminae, each of width Wm ( 5nn= WmN ) , 
the contact deformation at the axial mid point of a given laminum 
is r 


(Stcnrij ="Daij “ C(n Wm(k " '2) "Vm 


-CD 


Km 


7^1 = 0 tJj,= I 


3. 1. 1.5 


^ 4 and 2 are angular displacements (misalignment) of the 
bearing inner raceway relative to the outer raceway. Figure 5 
illustrates the components of Eq. (^3.1. 1.5.), The variables iTm^ 
are the relative normal displacements of the roller axis toward 
the bearing raceways. 

At a given azimuth, , contact deformation varies with 
laminum location in the contact zone. The interference at the 
midpoint of the contact is a function of bearing inner ring 
displacements relative to the outer ring and may be equated 
to the sum of the inner and outer contact deformations. 



3. 1.1. 1.2 EHD Lubrication 


3. 1. 1.6 


The lubricant is assumed to be a Newtonian fluid which 
implies that the shear stress, viscosity, and velocity gradient 
are related as follows; 





3. 1. 1.7 


Since, in EHD contacts, the lubricating film is extremely thin, 
the lubricant velocity gradient is assumed linear. Therefore, 
Eq. 3. 1.1.7 becomes: 


'V 


3. 1. 1.8 
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Where v is the roller/raceway or roller/flange sliding velocity, 
and h is the film thickness. For most lubricants the viscosity- 
temperature-pressure relationship can be presented in the form 
of Figure 6 . For a given temperature, an approximation to the 
curves similar to those of Figure 6 is: 

^ ^ 3. 1 . 1.9 

Where O' represents pressure and Cj, C 2 and C 3 are constants 
calculated from input pressure-viscosity data using a least 
squares fit procedure. 

For a given laminura, the pressure distribution in the y 
direction is given by: 

With y measured from the center of the contact. 


t = y/b 




3.1.1.12 




2 




r refers to the roller and raceway crown radii. 


3.1.1.13 


At the roller-end-flange contact, (REFC), the contact shape 
is assumed elliptical and the pressure distribution is assumed 
to form an ellipsoid. At any point within the contact pressure can 
be given by: 
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^ ZTTo^Uj 


3.1,1.14 


Where : 


Si- 


3.1.1.15.1 


t - Vw, 


3.1.1.15.2 


Given the principal radii of the flange and roller end surfaces; 


for the roller; 




3.1.1.16 


for the flange; 


3.1.1.17.1 





3.1.1.17.2 


the semi-major ( a j ) and semi-minor (bj) axes of the contact ellipse 
can be determined according to the methods of ref, 3. Using 
Eq. (3. 1.1.9, 10 and 14) viscosity can be expressed as a function 
of pressure. 


The sliding velocity, T/tm] at any laminum in a raceway 
contact, depends upon roller speed, raceway speed and geometry 
of the deformed bodies. Palmgren, ref. 10 states that 2/3 
of the contact deformation occurs in the roller and 1/3 in the raceway 
body. Therefore, 


n/ir ■ - -L 


- (D<- ^ Wjl 3.1.1.18 


The sliding velocity at the center of the REFC is assumed 
to act at each point within the contact ellipse and is given by: 


ATfj. - (a -H'jioj 
11 


3.1.1.19 
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and H are defined in Figure 2 


~ ^c. 


3.1.1.20 


3. 1.1. 1.2.1 EHD.Filtti Thickness 

EHD films are assumed to be developed in the line contacts 
at both raceways and at the point contact of roller end flange. 

3 . 1 . 1 . 1 . 2 . 1 . 1 Raceway Film Thickness 


According to Dowson et al (11), the lubricant film thickness 
in a rolling line contact is defined by: 


o.U 



3.1.1.21 





3.1.1.22 


3.1.1.23 


3.1.1.24 


= D/e 3.1.1.25 

and Xo are based on lubricant bulk temperature. 

Others have developed formulae similar to the foregoing; however, 
Eq. ( 3 . 1 . 1 . 21 ) is used herein. 


3. 1 . 1 . 1 . 2 . 1 . 2 Flange-Roller End Film Thickness 


For an elliptical point contact Archard et al ref. (12) have 
developed Eq. (3 . 1 . 1 , 26) to describe the lubricant film thickness 
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X. 


?■?... U . -Fj 

.1 1 . 






.a 4= 


EL 

0.1^ 






tSJ 


3.1.1.26 


= Vjj. 4- (p/7.- m) 


3.1.1.27 




3.1.1.28 


+ S.g^T 

^ tsPH 




can be 


For use in Eq. (3 . 1 . 1 . 21 , 22 and 26), the terras 
e developed from Eq. (3 . 1 . 1 . 9). ' 


3.1.1.29 
and X o 




^:0 


r ^ 
= c, e. 




Ac 


-M = 


= c>C^e 


-\ C.-L 


3.1.1.30 


3. 1. 1. 31 


2 \! Cx 

3. 1 . 1 . 1 . 2. 1 . 3 Raceway Traction 


Using Eq. (3. 1.1.8, 9, 10, 18 and 21), the raceway traction 
forces are obtained by numerically integrating Eq. (3.1.1.32). 




Fr^j - t, 

I Ji 


3.1.1.32 


tc--l 


3. 1.1. 1.2.4 Roller End Flange Traction 


Using Eq. (3. 1.1.8, 9, 14, 19 and 26) the roller end flange 
traction force is obtained by numerically integrating Eq. (3.1.1.33) 
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'tfj 


3.1.1.33 


F^j is assumed to act at the center of the contact in the direction 
of sliding at that point. 


3. 1.1. 1.3 Lubricant Hydrodynamic Drag Force 

In fluid-lubricated roller bearings, the rollers orbit the 
bearing axis, moving through space filled with a lubricant-gas 
mixture. This mixture and orbital motion effect a drag force on 
the roller which can be described empirically as follows: 


Foj = f — ^ 3.1.1.34 

Approximate values of the drag coefficient can be obtained from 
ref. (13). In Eq. (3.1.1.34), fluid density p is the mixture 
density. It is approximated by the quantity (lubricant weight 
contained within the bearing boundary dimensions) -f- (the free 
space volume within the boundary dimensions), expressed as a 
decimal fraction, designated P. 

3. 1.1. 1.4 Roller-Cage Normal Force 

A roller exerts a force Otj on the cage, if in the absence 
of the cage web constraint, it would tend to orbit the bearing 
axis faster than cage speed, CJg. Conversely, the cage pushes the 
roller with a force if in the absence of the cage web constraint, 
the roller orbiting speed tends to be less than cage speed. As 
noted earlier, the roller-cage forces are neglected as are any 
roller-cage web forces which would result from skew moments. 

3. 1.1. 1.5 Inertial Forces and Moments 

Eq. (3.1.1.35) is used to calculate roller centrifugal force 
at the roller center of mass. 


G Fj - 2^ rn 0. 


3.1.1.35 
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m is the roller mass. 

Roller gyroscopic moment in the XZ plane is given by: 



r 


3.1.1.36 


where J is the roller mass moment of Inertia. The gyroscopic 
moment is counteracted by the roller raceway contacts. 


3, 1.1. 1.6 Skewing (Yaw) Moments 


Owing to axial variation of frictional stresses acting on the 
roller/raceway contacts and REFC friction, a skewing moment 
occurs about the roller centroid, tending to cause yaw or skew 
of the roller axis. This skewing moment is resisted by the cage 
web elements, the REFC normal load and the restoring moment 
generated at the outer raceway, or all elements in consort. 

The magnitude of the skewing moment is: 





3.1.1.37 


ms-Z • 

By allowing the roller to adopt a non zero skew angle, Eq. (3.1.1.37) 
may be equilibrated to zero. However allowing the skew angle to 
be non zero makes the roller equilibrium problem twice as complex 
as it is presently formulated. Since the skew moment is generated 
by friction forces which for well lubricated contacts are at least 
two orders of magnitude lower than their respective contact loads, 
exclusion of the skew moment from the roller equilibrium solution 
should not have a major effect on the accuracy of the total bearing 
kinetic solution. 


3. 1.1. 1.7 Equilibrium of Forces and Moments 

In «ach principal direction, the forces acting on a roller 
constrained by raceways, guide flange or cage equilibrate to 
zero. If the roller does not contact the cage web, an acceleration 
or deceleration results in the y-direction. From Figure 2 the 
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three equilibrium equations are: 


^ i Qffi I 


ro-i 

mo-Z 


Coi>°< ir\ ~ CFj > Cos°c^ = O 


Z 


^I0o<ni “ - O 

m^Z 


A^. 


04 


dif 


3.1.1.38 


3.1.1.39 


3.1.1.40 


When cage/roller contact is absent, Qe.\~0 ; when such contact 
occurs ^ojs-^c.- 

For fluid-lubricated bearings, friction moments about the 
roller axis are obtained by integration of the product of friction 
stress and the corresponding roller radius: 


K-tJ 


Hmj ' \:p«.oij Ji 


3.1.1.41 


K.-\ 


Using Eq. (3.1.1.41) the following moment balance about the roller 
axis is obtained: 


m^Z 




m-i 




3.1.1.42 


3.1.1.43 


3.1.1.44 


Additionally, the moments in the XI plane must balance. 
Therefore : 


rn^ )f=>4 

y -(L-Lcrn) Si A ^ 3. 1.1.4. 
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Eq. (3. 1.1.6, 38, 40, 42, 45) can be solved for the unknown 
values of and 0 &] , if values of bearing ring 

relative displacements and cage speed are- known or assumed. 

3. 1.1. 2 Bearing Ring Loading and Displacements 

3.1. 1.2.1 Applied Loading 

As shown by Figure 1 axial (thrust) loading F j can be, 
applied to the bearing along the X-axis. Additionally, radial 
load applied in the YZ-plane can be resolved into components 
F2 and F 3 acting along the Y-axis and Z-axis respectively. If 
axial load is applied at a distance from the X-axis, an over- 
turning moment results, the components of which are F^ and F^ 
in the XZ plane and XY plane respectively. Similarly, moments 
resulting from eccentrically applied radial loading can be 
resolved into F^ and Fg components. 

3.1. 1.2.2 Displacements 

Owing to applied loading, the bearing inner ring is displaced 
relative to the outer ring as shown by Figure 7. Five displacements 
Al corresponding to the five loading directions are required to 
define the position, assumed by the inner ring, necessary to 
achieve an internal distribution of load which equilibrates the 
applied loading. 

Interference 


Figure 5 illustrates roller bearing interference at a given 
azimuth . 


Az] = Ai Ai Co 




3 . 1 . 1.48 
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is the axial interference, Az.^ is the radial interference 
and ky is the total ring roller interference at the center 1 am in urn 
according to Eq. (3.1.1. 6^ 

3. 1.1. 2. 3 Equilibrium 

The bearing applied loads, Fi, must equilibrate the summation 
of the individual rolling element loads according to Eq. (3.1.1.49 
53). 


J-Z 



3. 1. 1.49 

_Fz*- ^ (Oz-i Co^‘<2.- 0^4 

3.1.1.50 

F*, - 2 CQzJ - Q(j Co-s 

3.1.1.51 

4 -B 

-1 

J-2 

[e 2 - 

Viz. Co 

L Co-t Qi;.^ Co *b 

3.1.1.52 

4=1 


3.1.1.53 


+ - Qi'lsvnfi 

3. 1.1.3 Cage Speed 

3.1. 1.3.1 Cage Rail/Rinq Land Force 


Four basic cage types are used in the tapered roller bearing: 

1. Roller riding (RR) 

2. Inner ring land riding (IRL) 

3. Outer ring land riding (ORL) 

4. Inner and outer ring land riding "Z" configuration. 
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These are illustrated schematically by Figure 8. RR cages are 
usually not used in critical high speed applications. The choice 
of an IRL, ORL or "Z" cage depends largely upon the application 
and designer preference. With inner ring rotation, an IRL cage 
is driven by a force between the cage rail and inner ring land as 
well as by the rolling elements. ORL cage speed is retarded by 
cage rail/outer ring land drag force. The "Z" cage is driven at 
the inner ring and retareded at the outer ring. The magnitude 
of the drag or drive force between the cage rail and ring land 
depends upon the resultant of the cage/rolling element 
loading, the eccentricity of the cage axis of rotation and the 
speed of the cage relative to the ring on which it is piloted. 

If the cage 'rail/ring land normal force is substantial, hydrodynami 
short bearing theory (Ref. (14)) might be used to establish the 
friction force F^j^. For a properly balanced cage and a very small 
resultant cage/rolling element load, which are assumed here, 
Petroff's law can be applied to both cage^ rail-land surf aces ; e, g. , 

FcUfj = 1^ 0 (^Jc~S2m ) 

1 - 3.1.1.54 

Ot 

where D 2 is the larger of the cage rail and ring land diameters 
and is the smaller. 


Fcl = E ^■cun 3.1.1.55 

n-i 

n = 1 refers to the roller large end side 
n = 2 refers to the small end side 

3.1. 1.3.2 Cage Torque Balance 

For steady state operation, the sum of the torques acting on 
the cage is ‘zero. Therefore: 

Qcj Feu — O 3.1.1.56 

y-' 
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3. 1.1.4 Bearing System Solution 

For tapered roller bearings, Eq. (3.1.1.49-53) and Eq. (3.1.1.56) 
can be solved for the displacements tk[, and cage speed, 6)(. • 

This solution is obtained for a given set of bearing and lubricant 
temperatures and bearing diametral clearance or end play. As 
demonstrated in section 3.1.3 diametral clearance is a function 
of bearing operating temperatures, ring speeds and rolling 
element loads. It is essential that (1) the temperatures 
assumed for operating clearance, lubricant film thickness and 
traction calculations be consistent with actual operating temperatures 
and (2) that assumed clearance based on assumed rolling element 
loads be consistent with actual internal loading. 

If calculated temperatures and rolling element loads are not 
consistent with the original assumptions, then the entire calculation 
or procedure must be repeated using new values of temperature based 
on the last set of calculated temperatures. This process is 
continued until a consistent set of kinetic and thermal conditions 
has been determined. 

3.1.2 Temperature Mapping Heat Dissipation Model 

At high operating speeds the performance characteristics 
of rolling element bearings tend to deviate substantially from 
the characteristics which would be expected based on an extra- 
polation of slow speed behavior. Cage and rolling element 
speeds tend to be lower and bearing frictional heat generation 
tends to be higher than predictions based on slow speed extra- 
polations. These effects are, at least in part, caused by 
hydrodynamic forces generated as the lubricant is pushed out of 
the roller paths and/or is forced to move with the roller set at 
high speed. These lubricant related forces are highly viscosity 
dependent and viscosity in turn is highly temperature dependent. 

The temperature dependence of viscosity is important in both 
hydrodynamic (HU) trad elastohydrodynamic (EHD) lubrication. 

Accurate temperature predictions result in more accurate 
viscostiy predictions and thus more accurate friction force and 
bearing performance predictions. 

Temperatures have an additional effect on bearing performance. 
Radial temperature gradients across the bearing tend to alter the 
bearing diametral clearance and thus the distribution of loading 
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among the rolling elements. Axial temperature gradients can cause 
a pinching of the roller at one end, and thus alter the distribution 
of load over a roller profile. 

The effects of bearing assembly and lubricant temperatures 
should be an integral part of a bearing dynamic performance analysis 
when temperature gradients are large. The tapered roller bearing 
dynamic analysis has provision to account for both of these 
temperature effects. The operating geometry analysis Section 
3.1.3, accounts for the effect of temperature on component 
dimensions and lubricant viscosity is determined for local 
lubricant temperatures. 

Most bearing analysis precedes bearing manufacture and most 
bearing tests are not instrumented sufficiently to yield all of the 
temperatures required for accurate analysis. To allow the 
extrapolation of experimental temperature data or to allow 
temperature estimates to be made without experimental information, 
a temperature mapping heat dissipation analysis has been developed. 
The analysis uses numerical methods of finite differences as 
demonstrated by Dusinberre ref. (15). System temperatures are 
predicted using the frictional heat generation rates calculated by 
bearing kinetic analysis and consider the following system 
variables ; 

• Bearing assembly dimensions and materials 

• Ambient air temperature and velocity 

• Circulating fluid coolant through the housing 

• Lubrication scheme 

The numerical method of solution requires the bearing system 
to be divided into discrete mass segments, each represented by a node 
point. A map of the nodal points is presented in Figure 9. Figure 
9 represents a plane through the bearing axis of rotation, and it 
is assumed that all the system temperatures may be represented 
by such a single plane, i.e., the system is symmetric about the 
bearingaxis. High bearing speeds and large lubricant mass flow 
rates make this assumption reasonably accurate. This assumption 
implies that each temperature node is a circumferential ring 
around which no heat may flow. 
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3. 1.2.1 Lubricant Distribution Analysis 

The lubricant distribution model is not purported to be 
the result of a rigorous mass transport analysis. It is a 
starting point which Can be modified and improved as more 
definitive data become available. It is based on simplifying 
assumptions, many of which are intuitive. Although the distribution 
model was developed primarily to give better definition to the 
thermal analysis and provide a method to optimize a lubrication 
scheme, it was believed that the concept would prove useful in 
additional areas of bearing analysis including starvation and 
transient analyses in which the lubricant supply is terminated. 

In high speed, fluid lubricated, rolling element bearings 
practically full EHD films exist and all frictional heat generation 
is assumed to result from EHD and HD lubricant shearing. To 
effectively model the bearing heat dissipation and bearing 
temperatures by the finite difference method, the temperature 
grid must be sufficiently fine to represent the real situation. 

It is obvious that at least three lubricant temperature nodes 
are required to coarsely represent the heat generated by a tapered 
roller bearing, ome at each raceway and one at the flange. It 
was conceived that the lubricant might enter the bearing cavity 
at four distinct locations; through the inner raceway at both 
ends of the roller path and by jet to one or both sides of the 
bearing. Continuity requires at least one node to which all 
lubricants can flow for exit. Thus, eight lubricant nodes 
are required as a minimum to define a relatively coarse system. 

With more than two lubricant nodes which simply define entrance 
and exit, a model must be developed which apportions the lubricant 
distribution such that continuity is maintained. 

As finally conceived, the tapered roller bearing lubricant 
node grid has a maximum of 23 nodes, 6 within the bearing which 
act as the originating points for all of the generated frictional 
heat. Referring to Figure 9 and 10 these six nodes are numbered 
3, 4, 5, 12, 13 and 14. Three nodes represent lubricant entrance 
to the bearing cavity via jets, nodes 8 and 11, and via the scoop, 
node 20. The scoop delivers lubricant through flow paths in the 
shaft and inner ring to both sides of the inner raceway at nodes 
15 and .17. Nodes 9 and 10 represent the lubricant which is jetted 
to but does not enter the bearing. Heat can be convected from the 
roller ends to these nodes. Nodes 2 and 6 accept all lubricant 
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which is rejected by or exits from the bearing cavity. Nodes 
1 and 7 represent scavenges at which no heat convection takes 
place. 

By its nature the model is complex. In an effort to avoid 
further complexity, it was decided that the lubricant flow rate 
to each node would be defined by specifying linear combinations 
of the lubricant input rates. If this were not done but instead 
the rollers were allowed to transport lubricant from the outer to 
the inner race, an iterative solution would be required along with 
inherent additional complexity. It is necessary therefore that a por- 
tion of the lubricant at each and every input node, ^nodes 8, 11, 

15 and 17, be distributed to each of the inner raceway lubricant 
nodes 12, 13 and 14. All six bearing lubricant nodes will 
receive some flow since it is assumed that centrifugal force will 
drive the lubricant at the inner race nodes to their counterparts 
at the outer races, nodes 3, 4 and 5 respectively. The lubricant 
flow rate entering the bearing from each of the nodes 8, 11, 15 
and 17 is divided into three auxiliary components designated 
Pip and Me, for the scoop and two component s Pl<,p and & , for the jets. 

In a high speed bearing the major portion of, lubricant flow is com- 
prised of Mb . From the jet nodes 8 and 11 Mg, is defined to be 
the flow rate of the oil distributed to the nodes at the inner 
race as a function of the lubricant velocity leaving the jet. The 
jet velocity reflects the lubricant’s ability to penetrate into 
the bearing. From the through race entrance nodes, 15 and 17, 95 
percent of Mb enters the lubricant nodes at the edges of the inner 
race, nodes 12 and 14 respectively. Node 12 accepts all of the 
frictional heat generated at the roller end flange contacts. 

5 percent of Me, is assumed to exit the bearing. 

From node 15 and 17 tip is distributed to all six of the 
lubricant nodes within the bearing. From nodes 8, 11, 15 and 
17 n«,p is distributed to the lubricant nodes at the inner race, 
nodes 12, 13 and 14. Pl^p decreases with increasing cage speed. 

The effect of high cage speed is to increase the relative amount 
of lubricant distributed to the inner race lubricant node nearest 
the jet and results in corresponding decreases in the amount of 
lubricant reaching the more distant nodes. The same effect is realized 
with through race lubrication. In addition high speed results in 
less lubricant distributed to the inner race nodes in favor of 
the outer race. The lubricant distributed by the terms Mspand 
rip represent the lubricant transported by random splash and surface 
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wetting. These terms insure that at least a small portion of 
lubricant is distributed to each and every bearing internal 
lubricant node from each lubricant inlet source. 

, A detailed description of the determination of and 

Hf as well as the flow rates from node to node is presented below. 
Figures 10.2 and 10.3 give a graphical summary of the lubricant 
distribution assumptions for the lubricant supplied by under 
race cooling and for the lubricant supplied by direct jetting 
respectively. 

3. 1.2. 1.1 Under Race Lubrication 

As indicated in Figures 9 and 10, the under race lubrication 
is admitted to the bearing via a scoop and a set of flow paths; 
the number of channels which comprise given flow path is 
specified at input. The scoop efficiency, i.e., that fraction of 
the lubricant input to the scoop which actually enters the 
bearing, must be input to the analysis. When both under race 
flow paths are used, the ratio of the lubricant velocities flowing 
through paths 1 and 2 is obtained by equating the ratios of the 
pressure heads developed by the centrifugal pump effect, to the 
ratios of the negative pressure heads attributable to friction: 


V. ^ T - "Rp-iOTeAisiCE ) 

Where V refers to the lubricant velocity 


3. 1. 2. 1 


D refers to the flow path channel diameter 
R refers to the radial distance of the flow path entrance 
and exit points from the shaft centerline. 

L is the total length of each respective flow path. 


Knowing the ratio of the velocities at which the lubricant is 
flowing, the ratio of the lubricant mass flow rates may be obtained, 
referring to Figure 9, 10.1 and 10.2. 




3. 1.2.2 
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The designation refers to the lubricant mass flow rate from 

node i to node j. and l^l- 2 .e,xs- may be obtained from 

Eq. (3. 1.2. 2 and 3). N refers to the number of radial channels 
spaced around the circumference for each flow path 


H 7.0/ lA 

M s i s the total mass flow rate to the scoop. 


3. 1.2. 3 

3. 1 . 2. 4 

From Figure 10 


^0 Hs 3. 1.2.5 

HiOjfc is the rate at which the lubricant is rejected from the scoop. 
Continuity requires, 


'S’ - 1^7.0, 7-S 


3. 1.2.6 


and 



3. 1.2. 7 


From nodes 15 and 17, the under race lubricant is assumed to be 
distributed symmetrically to the bearing surfaces such that; 


ni€,iz- __ 
1 



M 77j 





NviO 


\ 

HHj 14 


77j 'S 
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Therefore, is fully defined if Plis^j is fully defined. 

K and j are symmetric about a surface through nodes 4 and 13 . 


What follows is an explanation of all nic,^ mass flow 
rate terms. To ease the explanation, we will temporarily define 
a term M equal to M 22 ,i 5 a speed dependent variable S: 




4 - 


e. 4 - X 3 


3.1.2.10 


The purpose of the variable S is to account for the effect 
of speed and the centrifugal force field. High speed tends to 
limit the amount of lubricant which is transported from one side 
of the inner race to the other. Due to centrifugal force, the 
major portion of the lubricant is flung to the inside surface of 
the outer ring and from there flows out of the bearing cavity. 

As noted earlier, we will define several auxiliary terms: 


■a 

Hp-5. 0.-Z7 H (,l-s) 


3.1.2.11 

3.1.2. 12 


3.1.2 . 13 


The lubricant flow rates from node 15 to the other nodes 
within the system are assumed to be given by the following 
linear combinations of the quantities nsp, andna with the 
nodes nearest node 15 receiving proportionately more lubricant 
than the more distant nodes. Referring to Figure 10 ; 




= 0 . 0 s Ha. 

CM 

• 

• 

CO 

1 

I 


Hp 

3.1.2.15 

To the outer raceway J 

H 

I 

Hp 

3.1.2.16 

1 

L tSj s 

- 0 . U 1 

3.1.2. 17 
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63^1 Hf +o.s'H"i.p -t o-^s Me, 


To the inner raceway 


= 0.(12- ■‘rd. S2»5 

« t I 

Ki>;(4 = (O.ass o.\(o1 K-sp 


3.1.2. 18 


3.1.2. 19 
3.1.2.20 


Examination of Eq. ( 3 . 1 . 2 . 10-20) and Figure 10.1 and 10.2 
reveals that at high speed, the major portion of the under race 
lubricant entering the bearing cavity at node 15 flows to the large 
end flange, node 12, and then to the outer race. Relatively 
small portions of the lubricant are assumed to be splashed or 
transported to the inner ring surface. As noted earlier, the 
lubricant flows from nodes 15 and 17 are assumed symmetric one 
with the other, thus Eq. (3.1.2.10-20) also describe the under 
race lubricant distribution from the roller small end side of the 
inner race. 


3. 1.2. 1.2 Jet Lubricant Distribution 


Just as the lubricant distribution from nodes 15 and 17 was 
assumed symmetrical, the distribution from the jets, nodes 8 and 
11, is also assumed symmetric. To reduce repetition, only the 
equations for the distribution from node 11 at the roller small 
end side of the bearing will be developed. 

The lubricant is assumed to be jetted between the cage and 
the inner ring. However, the lubricant's free access to the bearing 
cavity is limited by two considerations; 

1. The radial gap between the cage and inner ring may be 
smaller than the jet diameter, and 

2. A portion of the annular gap between the cage and inner 
ring is occupied by the roller ends. 

In Eq. (3. 1 . 2 . 2p) the first and second terms account for these 
respective considerations. 
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It is assumed that the jet is directed at a point midway between 
the inner raceway and the pitch diameter, is the jet efficiency, 

defined as the ratio of the amount of lubricant which actually 
enters the bearing, po the amount which is jetted to the bearing. 

Dj refers to the jet diameter. is the radial gap between the 

cage rail I.D. and the inner ring and O.D. Jet efficiency can 
also be specified at input, in which case Eq. (3.1.2.20)is not 
used. 


The set of auxiliary variables is used to describe the 
distribution of the under race lubricant and will be redefined in 
terms of the jet lubricant flow rate M which enters the bearing, 
(Refer to Figure 10.3). 3 1 2 22 


HspH- ^o.\e>A-o.szS)ri 
where S is defined in Eq. (3.1.2.10) 


3.1.2.23 


Ha = O.'S-Z. (,|-s)M 3.1.2.24 

% 

The lubricant flow rate Me. is distributed across the 
bearing inner race according to a function of several variables. 

1. The lubricant velocity leaving the jet 

2. The cage speed and cage type 

3. The bearing geometry 

The time available for the jet to penetrate between rollers 
is given by 


© =» 


r 


•Z.TT 




0.435 O 
e - D/-Z. 


3.1.2.25 


A secand variable is a ratio of the distance the jet can travel 
during the time ©, to the length of the roller r’ • 


q? ^ Vj S> 3.1.2.26 

L 
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where Vjis the velocity of the jet; 

1 / - . 

3.1.2.27 

P is the lubricant weight density 

is the cross sectional area of a jet at the roller small end side 
K 3 is the number of jets at the roller small end side. 

If the cage is outer race guided, 'f* is limited to 1. If the 
cage is inner race guided,^ is limited to 0.666. The distribution 
of the lubricant from the jets at the roller small end side of the 
bearing to the nodes at the inner race may now be determined as 
follows : 

For 0.75 < ^ <1.0 

3.1.2.28.1 

3.1.2.29.1 

3.1.2.30. 1 


= 0 , 

= (( - he> 

for 0.0 < < 0.25 

I ^ 

Miijii. - o 

= Me 
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3.1.2.29.3 
3. 1.2. 30.3 


3.1.2.28.2 

3.1.2.29.2 

3.1.2.30.2 


H ^ ^ i ) H 
M 11/?) *1a> 

M n/4 - ^ ha 
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The flow rates in the three equations Eq. (3.1.2.28, 29 and 
30) are supplemented by additional flow rates which are not 
related to the velocity of the lubricant leaving the jet. This 
is done to insure that some lubricant reaches nodes 12 and 13 
even if the Jet velocity is insufficient to penetrate to these^ 
nodes. These flow rates are assumed as follows: 


and 







VI 






3.1.2.31 

3.1.2.32 

3.1.2.33 


3.1.2.34 


The lubricant flow rate which is j.etted to, but does not 
directly enter the bearing, is 


H„,,. 

3. 1.2. 1.3 Free Lubricant Distribution 


3.1.2.34 


The distribution of the lubricant from the jets and from the 
under race entrance points has been examined. It remains to 
define the flow rates from the inner ring to the outer ring nodes 
and finally from the outer ring nodes to the scavenges. 


All of the lubricant which has flowed to a given inner race 
node is assumed to leave the node by flowing radially outward to 
the corresponding node at the outer race. 


RESEARCH 






j 
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The flow of the lubricant out of the bearing along the outer 
race is governed first by the presence or absence of scavenges. 
Referring to Figure 10, if scavenges exist at both sides of the 
bearing, the lubricant flow direction for a tapered roller 
bearing is assumed, from node 5 to 4 to 3 to 2 to 1 and from 
6 to 7. If on the other hand, scavenges exist only at one side 
of the bearing, all of the flow is toward that side. The 
scavenges are represented in the system by nodes 1 and 7. 

It is significant to note that the rate at which lubricant 
enters a given lubricant node can be determined simply by: 


2 3.1.2.39 

j -I 

All Mjt terms not specifically defines are set to zero. 

3. 1.2. 2 Heat Dissipation and Temperature Mapping 

At steady state conditions, heat generation rates are constant 
and for a given temperature node, heat flow to a node must equal 
heat flow from the node. Referring to Figure 9; 




t. = \ 

•L‘6 -Is 7. 


3.1.2.40 


QjjC. refers to heat flow from neighboring temperature nodes. 
The subscripts j,L indicated flow from node j to node i. 

Ql refers to the frictional heat which is generated at node i 
refers to heat transfer with ambient air. 

t , 

All Qjt and terms are temperature dependent. This 

enables the system of equations, developed from Eq. (3.1.2.40) to be 
solved for the individual nodal temperatures. 

3. 1.2. 2.1 Heat Sources 


As noted earlier, the heat sources in this model are all 
assumed to result from hydrodynamic or elastohydrodynamic friction 
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forces, i.e., lubricant shear. 

A segment of the bearing containing one roller contains five 
sources of frictional heat generation. Three of these are calculated 
on the basis of EHD lubricant shear. The hydrodynamic friction 
results from the rollers plowing through the lubricant and also from 
the cage rotating centrically on the inner and/or outer ring land. 

3. 1.2.2. 1.1 EHD Frictional Heat Generation Rate 

For a given roller in a tapered roller bearing there exists 
three EHD lubricated contacts. The frictional heat generation 
rates in these contacts may be obtained by combining Eq. (3.1.1.18 
and 32) for the raceway contacts: 






(ts-Kl 


vn 


2 b 








3.1.2.41 


' - w 

= 1 refers to the outer raceway contact IT) = 2 to the inner 
raceway. Eg. (3.1.1.19 and 33) may be combined to obtain the 
heat generation rates for each roller end flange contact. 

4, 


= q<.iW.iVyj f I di j 


J. J 


3. 1.2.42 


In the above equations J converts heat generation rates in 
Ibf-in/sec to BTU/Hr. 

The assumed axial symmetry allows heat to be input to the 
model as though the heat were generated at a single segment, 
such that for the three EHD contacts; 


A"*-) 




fin = 1 Outer Race 
tn = 2 Inner Race 
= 3 Flange 


3.1.2.43 


3. 1.2. 2. 1.2 Roller Drag 


The. rollers, as tbey orbit within the bearing cavity, displace 
and shear the lubricant and thus generate heat according to the 
semi-empirical relationship from Ref. 4 (refer to Eq. 3.1.1.34). 
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3.1.2.44 

If there is a substantial amount of lubricant within the 
bearing cavity, and the cage orbital speed is high has been 
found to be several times higher than the sum of all heat 
generated in the EHD contacts. 


3.1. 2. 2. 1.3 Cage Drag 


Assuming negl igible. cage rail to ring land force, the 
frictional heat generated at the cage rail-ring land contacts 
may be estimated by Petroff's equations as follows: 


A fjo C I ~ i ^^0 A el n 

n-\ 


3.1.2.45 


where the subscript 1 refers to the roller large end side of the 
bearing and 2 refers to the small end side. A is the area at 
the cage rail/ring land interface. d is the land diameter and 
^d is the radial gap between the cage rail and ring land. 


3. 1.2. 2. 2 Frictional Heat Distribution 

The lubricant churning heat Eq. (3.1.2.44) is assumed to be 
generated primarily as the roller set pushes the bulk lubricant 
around the raceway, thus the churning heat is distributed 

to the three lubricant nodes at the outer race. Also, it is 
assumed that the amount of churning heat distributed to each of 
the three lubricant nodes is proportional to the amount of 
lubricant at each of the three nodes. The total lubricant per 
hour at a given outer raceway node j is given by: 


t 


Hi ^ 2 VI M 


3.1.2.46 


A proportionality constant Kn can be developed 


K.n = 


Mn 


n “ -s^ 


3.1.2.47 


Z rt. 
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Heat from all sources is distributed according to Eq. (3.1.2 
52.2) 


Qa ~ 1 ^ V ^4 3.1 

0 Wx. 3 . 1 

For an outer ring land riding cage. Fig. 8; 

.. .. -Qit =■ U 4 ■♦''4 (S^ 1 + Z "SjS 3.1 

= Ws- ^ 2- 3.1 

3.1 

For an inner ring land riding cage: 

= 'A R + ICk 3.1 

jQ \o- = ^ '5 4 W x-V U 5 3.1 

Q W ■». 4 Vz t4 s 3.1 


For the "Z" configuration cage is defined as follows; 

G ?3 by Eq. (3.1.2.50.2) 

i * 

Q S by Eq. (3.1.2.50.1) 
by Eq. (3.1.2.51.2) 

Qh by Eq. (3.1.2.52.1) 
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3. 1.2. 2. 3 Heat Transfer 

3. 1.2. 2. 3.1 Ambient Convection and Radiation 


For the housing external surface, considering convection 
and radiation ref. (16) gives, 


+ p.5‘‘45 DvJ e. 




\ \(90 J V \oo 




3.1.2.53 


T - refers to temperature in degrees Farenheit 
w - refers to the nodal cylindrical surface width-feet 
D - refers to the diameter of the surface-feet 
e - is the housing emmissivity. 

- is the heat flow rate in BTU/Hr. 


Eq. (3.1.2.53) may be modified by miltiplication of the 
entire right hand side by ^ to account for forced convection 
at the housing. According to the data of ref. (17); 


^=■14- 


\/(V 

-r 
I tc 


3.1.2.54 


where is the air velocity in ft/sec. 


3. 1.2. 2. 3. 2 Conduction 


Between adjacent nodes in the (bearing, shaft, housing) 
assembly, heat is transferred by conduction such that; 



3.1.2.55 


is the area of normal to the heat transfer path between 
node s i and j, f t2 

Ll'j is the linear distance between nodes i and j, ft 
K is the thermal conductivity of the material between nodes 
i and j, BTU/ft Hr®F 
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If, between nodes, heat must flow serially through two 
distinct materials having thermal conductivities Ky and K 2 . and 
thicknesses Lj and Li then an equivalent heat transfer coefficient 
must be determined such that', 

3. 1.2,56 

Ai.; = Km ( \ 

V 3.1.2.57 

Atj 1^1 V<-T- 

Note: /t] ^ '^Tr\ and > and the distance between nodes 

i s L I ‘j - U 1 4 - L 7 _ 

If, between nodes, heat must flow through parallel paths of 
two distinct materials having thermal conductivities and 
and areas normal to heat flow of and A^, then an equivalent 

heat rransfer coefficient/^ is calculated such that; 


Xc) = (k. 




3.1.2.58 


Note; that the cross sectional 
3. 1.2. 2.3. 3 Convection 


area normal to, heat flow is 


ao 





Eq. (3.1.2.59) describes the convective heat transfer between 
two adjacent temperatures nodes, one in the assembly, node j, 
and One in the lubricant or fluid coolant, node i. 





3.1.2.59 


h - is the convective heat transfer coefficient. BTU/(hr 
ft2°F). 

Tm - is the mean temperature of the lubricant at node i 
defined as follows: 
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^ . lituTk 
Tl ^ 






Z Ht^'u 


it -I 


3. 1.2.60 


A - is the heat transfer area normal to heat flow ft^. 


Three relationships for the convective heat transfer co- 
efficient are used in this analysis. 

1. Laminar flow in the coolant and in each under race flow 

path. 

2. Laminar flow in the housing and on the shaft. 

3. Turbulent flow within the bearing cavity. 

3 . 1 . 2 . 2 . 3 . 3 . 1 Laminar Flow in Tubes 


For Laminar flow within the helically wound outer ring coolant 
tube and also for the under race lubricant flow paths, Eq. 
(3.1.2.61) from ref. (18) is adapted to describe the heat transfer 
coefficient. 







3.1.2.61 


Where D is the flow path diameter, ft 

L is the flow path length, ft 
M is the fluid mass flow rate 
Cp is the fluid specific heat 

T) is the fluid dynamic viscosity. The subscript b refers 
to the evaluation of the viscosity at the fluid bulk tem- 
perature W refers to wall temperature. 

K is the fluid thermal conductivity BTU/Hr-ft®F 


Of the terms within Eq. (3.1.2.61) K, cp and f) are temperature 
dependent. However the temperature dependenceof K^and Cp is 
weak, therefore, ,the values of K and Cp are determined for the 
lubricant based only on the inlet temperature. From ref. (19); 


C.p 


0.1? T* 


3. 
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r 




o. o (eT{(o 


I - . poo -h. 


3.1.2.63 


Where i s the lubricant specific weight at 60°F. 

Eq. (3.1.2.62 and 63) are valid for typical lubricating oils. 


Since the housing coolant may be any fluid, the coolants 
thermal conductivity, specific heat, and viscosity are input 
to the model. 


To simplify Eq. (3 . 1 . 2. 6l), the term is set to 

1.0. This is justified since Tjj and T^^ will not be substantially 
different. Also, with such a small exponent on the viscosity 
ratio, even if Tj^ is several hundred degrees higher than Tj^ such 
that ) is 5, is only 1.25. 

3. 1 . 2. 2. 3. 3. 2 Laminar Flow Unrestricted 


There are nodes within the assembly adjacent to which the 
lubricant is not forced to move at high speed. Referring to 
Figure (9), these nodes are numbered 1, 2, 6, 7, 8, 11 and 20. 
At these nodes, ^an adaptation of the convective heat transfer 


for laminar flow over 


h - Q.S'bT.KFy 


a flat plate is used such 

it 


coefficient 
that : 


3.1.2.64 

gn must have units of ft2 when K has units of BTU/Hr-ft°F 
Wp Vust be set to its absolute value. Pr represents the dimensions 
prandtl number, defined as: 



Fr 




must be expressed in units of Ib/hr-ft. 


3.1.2.65 
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3 . 1 . 2 . 2 . 3 . 3 ■ 4 Turbulent Flow 


Referring again to Figure 9 , the lubricant at nodes 3, 4, 

5, 12, 13 and 14 is forced to travel over the bearing internal 
surfaces at relatively high speeds. The bulk lubricant is assumed 
to be carried between rollers at the cage rotational speed. 

Determining the Reynolds number for a 120 mm bore tapered 
roller bearing lubricated with a MIL-L-23699 lubricant at 363°K 
(200°F), rotating at a shaft speed of 12500 rpm, a Reynolds number 
may be estimated from Eq. (3.1.2.66) to be, 7.5 X 10^, well into 
the range of turbulent flow. 

y *3 ^ 5 3.1.2.66 

In the turbulent region Kent, ref. (16), gives Eq. (3.1.2.67) to 
describe the convective heat transfer coefficient for oil. 


h = c 


V 


1 


• o.foi 


3.1.2.67 


V - is the relative lubricant/surface velocity ft/sec. 
has units of Ib/hr-ft 

C is a constant having a value of 122 if the lubricant is being 
heated or 91.5 if the lubricant is being cooled. 

h has units of BTU/h r-f t 2- ®F . Inspection of Eq. (3.1.2.67) indicates 
that typical convection heat transfer coefficients for oil in 
turbulent flow are of the order of 5000 BTU/hr-ft2-°F . Coefficients 
of this magnitude effectively help bring the bearing surface 
temperatures into near equilibrium with the adjacent lubricant 
temperature . 


Eq. (3.1.2.64) and Eq. (3.1.2.67) are viscosity dependent 
and as stated earlier, viscosity is temperature dependent. The 
temperature viscosity relationship assumed, is: 






3.1.2.68 
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T*^ it '' 

Where '(i and are input to the analysis. is the 
lubricant viscosity at temperature T^* and is the viscosity 
at T-i both®^tmo spheric pressure. ^ 


3. 1 . 2. 2. 3. 3. 5 Heat Flow by Mass Transfer 


Heat transfer between fluid temperature nodes occurs entirely 
through mass transport. It is necessary to develop an additional 
temperature , Ti , for each node in the lubricant, which is the 

temperature of the lubricant entering a given node. 

j-zs 


Til - 

2 Mji 

H • • 

Values of are determined by the lubricant distribution 

analysis. WithTxl. thus defined Eq. (3.1.2.60) may be simplified 


Tmi — ^ 

3.1.2.70 


Within the lubricant, the rate of heat transfer, from all nodes, 
to node i is given by such that 

(ijL ■=• Cp 3.1.2.71 

3. 1.2. 2. 4 System Boundary Conditions 

A portion of the system boundary conditions has already been 
mentioned i.e., the ambient air temperature and velocity. In 
addition, the temperatures at the shaft and housing ends must 
be specified. If the effective shaft and housing widths are suit- 
ably long, the end temperatures may be set to the ambient air 
temperature. If these widths are short, the end temperatures must 
be estimated. Specification of boundary temperatures in this 
way allows the bearing and lubrication system to be a heat sink. 
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The remaining boundary conditions are the lubricant and coolant 
mass flow rates and temperatures at input. 

3. 1.2. 2. 5 Heat Dissipation 

In addition to determining all assembly and lubricant 
temperatures, the thermal model also yields the rates of heat 
dissipation via: 

. Lubricant 

. Radiation and convection to ambient air 


, Outer ring housing coolant 
. To the shaft and housing end, nodes. 

3. 1.2.3 Significant System Temperatures 

For use within the bearing dynamic analysis four lubricant 
temperatures are required. 


The temperature of the lubricant at the flange is taken 
to be the temperature at node 12, at the outer raceway node 4, 
and at the inner raceway node 13, the bulk temperature, i.e. 
the mean lubricant temperature within the cavity is given by 
Eq. (3.1.2.73) 




_ 1 : M.T. 2 H.T, 


vt = a. 






3.1.2.73 


^ M. 


W,= ix 


As noted in Section 3.1.3 on operating geometry an effective 
temperature must be calculated for all assembly components. 

For instance, the equivalent temperature of the inner ring portion 
of section SI, refer to Fig. 9 and 11, is given by 












3.1.2.74 


<i. 


3. 1 . 2. 4 Closure 


The thermal analysis provides temperatures throughout the 
bearing assembly and lubricant. The analysis is easily adaptable 
to ball and cylindrical roller bearings. For all bearing types 
the analysis is meant to model high speed bearing operation 
where full EHD films are expected, the assumption of axial 
symmetry is valid. 
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3,1.3 Bearing Geometry Change Analysis From Cold Unmounted 
to Mounted Operating Conditions 

3.1.3. 1 Introduction 

As a bearing is taken from the shelf, mounted in a housing 
and on a shaft, turned up to speed and subjected to operating 
loads and temperature, the internal geometry of the bearing 
will change. To accurately analyze the dynamic performance of 
a bearing, its operating geometry must be known. An analysis 
has been developed to account for the following effects: 

1. Temperature changes and gradients. 

2. Initial and operating shaft and housing fits. 

3. Rotation induced ring radial growth. 

4. Uniform radial components of the rolling element- 
raceway normal loads. 

In conducting the analysis each bearing ring is conceptually 
divided into two parts by a radial plane perpendicular to the 
axis of rotation. These planes pass through the center of each 
raceway rolling path. By considering changes in roller path 
radii for the two sections the change in raceway contact 
angles and therefore the amount of roller pinching, as well 
as the change in bearing diametral clearance can be estimated. 

The basis for the major portion of this analysis is taken 
from Timoshenko, ref. 20 . The bearing rings are treated as 
thick walled circular cylinders of constant wall thickness 
subjected to the action of uniformly distributed internal and 
external pressure. 

The external pressure arises in the case of the outer ring 
from a press fit into the bearing housing. The internal 
pressure on the outer ring arises from the discrete roller loads 
which are regarded as uniform internal pressure acting on the 
outer ring. Similarly for the inner ring the press fit on the 
shaft provides a uniform internal pressure and the roller loads 
are regarded as a uniform external pressure. Fig. 11.1 shows 
the idealized sections used in the analysis. 
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3. 1.3. 2 Bearing Shaft and Housing Equivalent Sections 

Timoshenko's analysis requires that a ring be fully 
defined by specification of only internal and external radii. 

If two concentric rings are pressed together the analysis 
assumes that the rings are of equal width. Referring to 
Fig. 11.2 it is obvious that the outside surface of the 
actual inner ring section can not be defined by a single 
radius. However, a single radius may be found such that the 
cross sectional areas of the real and equivalent rings are 
equal . 

The equal area criterion may be applied to all four 
assembly aections SI, S2, HI and H2. Fig. 11.1 presents the 
equivalent sections for a tapered roller bearing assembly. 

Assembly refers to the bearing shaft and housing components. 

A second significant dimension associated with each of the 
four ring sections is the radius of the roller path at each 
section. Ultimately, it is the change in these four roller 
path radii, in addition to a change in the rolling element 
diameter, which permit a change in the bearing diametral clearance 
and a change in the inner and outer ring contact angles to be 
calculated. These radii are shown on Fig. 11.1 and are the 
radii to the actual inner and outer ring roller paths. 

The bearing geometry change problem requires the solution 
to similar problems for each of the four bearing sections. 

To help eliminate repetition and cumbersome subscripting, only 
section SI, will be covered thoroughly. However, differences 
in analysis which occur from section to section will be noted. 

3. 1.3. 3 Internal and External Equivalent Pressures 

Each of the four sections consists of either a solid or 
hollow shaft and inner ring and outer ring and housing respectively. 
The rings may or may not be in intimate contact. Intimate contact 
occurs when due to any of several factors, the outside radius 
of the inside ring is greater than the inside radius of the outside 
ring. A pressure is developed at that interface, which tends 
to expand the outside ring and collapse the inside ring. Although 
the pressure acting on the actual interfering surfaces are 
identical, the use of equivalent geometrical sections requires 
the use of equivalent pressure, such that the calculated pressure 
acting to collapse the inside will not be identical to the 
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pressure acting to expand the outside ring. 


The equivalent pressure concept allows the effect of unequal 
component width to be considered as well as the effect of shaft 
and housing shoulders. Both of these effects tend to increase 
the bearing seat stiffness above the value which^could be 
determined by considering only the shaft and housing diameters 
at the bearing seats. 


The effect of unequal inner ring and shaft width is the 
same for both sections Si and S2. The same is true for the 
outer ring and housing section HI and H2. This effect is 
taken into account by the factor X which applies to the pressure • 
where, referring to Fig. 11.2 

0,5- S (Zs - ^ 2, (3. 1.3.1) 

0.5 S (X. - $ 2. (3.1. 3.2) 

Where w denotes the element width and the subscripts I, 0, S 
and H refer to inner, outer, shaft and housing respectively. 


The effect on equivalent pressure resulting from shaft 
and housing shoulders is applicable at sections and housing 
shoulders is applicable at sections Si and H2 and is accounted 
for with the factor Y • Referring to Fig, 11,2 



^ 3 . 1 . 3 , 3 ^ 
(3. 1.3. 4) 



Where ^ is Poisson's ratio. 


(3. 1.3. 5) 
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Assuming multiplicative effects; 

Ym, - 7o„ 


( 3. 1.3,6) 


( 3, 1.3. 7) 


where denotes the pressure on the shaft OilD. as related to 
the pressure on the inner ring and To# 

similarly defined. The pressures (P) in Eq, (3. 1.3.6 and 7) 
are non zero only if an interference fits exists at the 
section in question. An interference fit may exist under cold 
mounted conditions but due to thermal gradients and high 
speed rotation, clearance may develop at operating conditions. 

To account for the elastic effect of the uniform radial 
components of the rolling element loads on the ring dimensions, 
an equivalent inner ring external pressure and outer ring 
internal pressure are calculated based on the minimum radial 
component of the rpllihg element loads, such that for the 
inner ring; 

-p _ 2 . (Qi.t C C?6 o<j, 

ztr d (3. 1.3.8) 

where Z is the total number of rolling elements 

is the smallest rolling element raceway 
radial load component at the inner raceway 

Wx is the inner ring total width. 

The external pressures thus calculated are assumed equal 
on both inner ring sections. An expression corresponding to 
Eq. (3. 1.3. 8) can be developed for the internal pressure 
acting to expand the outer ring. In fact, for a tapered roller 
bearing, will exist only when the applied load is predominantly 


(p. . COS i 
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axial or the bearing is preloaded, but the internal pressure 
acting on the outer ring, "Poi t will always exist to some 
extent, as a result of rolling element centrifugal force. 


3. 1.3. 4 Temperature Effects 

The effect of temperatures which occur at operating 
conditions is considered in two ways. A bulk effect based on the 
radii to the rolling paths (R) is calculated from Eq. (3. 1.3. 9) 

U ' ^ ^ ^3. 1.3.9) 

where 

U - is the change in rolling path radius 

^ is the component coefficient of thermal expansion, 

T - is the difference between the component effective 
temperature and 68°F. 

A similar expression may be written for the change in roller 
diameter. Referring to Fig. 9 for Section SI; 


^ (3.1.3.10) 


(3.1.3.11) 



The second effect of temperatures is the change in interference 
fit from cold to operating conditions such that the interference 
fit as a result of the initial fit and the temperature difference 
between the two components cai be calculated. Cold (293°K, 
(68®F)), shaft and housing fits Fcold are input to the 
analysis. A positive value reflects an interference fit. The 
analysis neglects asperity crushing, and, the fits for both 
sections of a given ring are assumed to be identical referring 
to Fig. 11 for section' SI; the fit (^hot) si. operating 

temperature is given by 


' F .. 




hot 


Si ' "T 5 ^ (3.1.3.12) 

is the fit considering operating temperatures. 
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3. 1.3. 5 Ring Radial Displacement Considering Surface Pressures 


At this point we consider the problem of evaluating thick 
walled cylinder, (ring) radial displacements as a function 
of the ring radial dimensions, physical properties, and uniform 
internal and external pressures. Timoshenko develops Eq, (3.1.3.13) 
for determining the ring/radial displacement (U), as a function 
of a general radius (r) to any point within the ring when 
con si dera ti on s of ring rotation are omitted. 

Vj * (3.1.3.13) 

— A 0. 


, _L c|y 
d7 r 3? 


The general solution to this equation is 


U k, f" 4- (3.1.3.14) 

r 

The expression for normal stress in the ring is given by: 





fk, (u t). k-l 



(3.1.3.15 ) 


The constants of integration Vi and can be determined 

from Eqs. (3.1.3.16 anH 17) which make use of the boundary 
conditions at the inner and outer surfaces of the ring. Note 
that the dimensions of inner ring sections Si, Fig. 11.1, are 
used in Eq. (3.1.3.16-19). These equations, however, are 
general expressions and are valid for all rings, 



"E p 

(3.1.3.16) 

<3.1.3.17) 
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Where E and ^ are the modulus of elasticity and Poisson's 
ratio respectively. The negative sign associated with Tx# 
reflect the sign convention, wherein a positive normal 
stress indicates tension. Solving Eqs. (3.1.3.16-17) for Ki 
and Kx and substituting the results into Eqs. (3.1.3.14 and 15) 
gives the general expressions for normal stress and radial 
deflection as functions of the internal and external pressures 


acting 

on the ring: 



(3.1.3.18) 

- b"*- 

(cx- b') 




. ' - j rb"*- 


-Pri 

(3.1.3.19) 


b " 




Using the appropriate dimensions and pressures, Eq. (3.1.3.19) 
can express the inward displacement of the shaft and the outward 
displacement of the inner ring when assembled with a interference 
f i t F 

HOT 


Ot b ~ 0 9b ~ ^ 


Hot 


(3,1.3.20> 


u 


Xb •= 




( 3.1.3.21) 


f - 5. 


J 


( 3.1.3.22) 


and are assumed zero. Using Eq, (3. 1.3.8) to obtain 

Tx, , the following expression may be developed forPcj which 
accounts for all variables in the interference fit problem 
exacept for the effect of ring and shaft rotation. 
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A- y \ / bc^Pi o\ 

fH„T + 1 I I C-- b-! ; 

H (^) * (§‘ ‘ ‘ “ 


Note, Pxi may not be negative, i.e. the effect of fitting one 
ring on another can not result in placing the common surfaces 
in tension. 

Having determined and ^To » Eq. (3;ii3.19) may be 

used to determine *Jx8 ( f , O by replacing the general radius (r) 
with the radius to the roller path, R. UxR (P.*) is defined 
as the change in rolling path radius resulting from: 


. The initial fit 


, The change in fit resulting from a temperature gradient 

. The effect of the rolling element load. 

3. 1.3. 6 Ring Rotation 

We must now examine the effect on the rolling path which 
results from high speed rotation. If the rotational speed is 
less than 100 rpm the rolling path radius is assumed to be 
unaffected. 

Timoskenko presents Eq. (3.1.3.24) to define ring dis- 
placement U in terms of a general radius (r), the weight density 
of the ring material, ( ^ ), and the ring angular velocity ( -^^ ) 

ci ^ ^ ^ ~ ^ ^3. 1.3. 24) 

Using the notation; 




pn: 


(3.1.3.25) 
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the general solution to Eq. (3.1.3.24) may be written: 


& >r 


U ---o-r 


( 3.1.3.26) 


The general expression for the normal stress is given by: 

^ ^ I" 5 /r-i 


E 


o 

o 



\ 




( 3.1.3.27 ) 


Using Eqs. (3il.3.26 and 27) and the principle of super- 
position, the effects of ring rotation can be considered. In 
this case, superposition allows a set of integration constants 
Ki and Kx to be calculated for a ring based on a change in 
pressure (P®**) at the internal and external surfaces. Eq. 
(3.1.3.27) may be written for the four specific surfaces of 
section SI, as follows; 


Inner Ring External Surface 
Inner Ring Internal Surface 


1 




_ P -Q (3.1.3.28) 


\ 


_ 


o ' 







= (3.1.3.29) 


Shaft External Surface 




l- 



U"*- 


J 


- -P^o (3.1.3.30) 


Shaft Internal Surface 






K). 


o 

e» 




2 5 0 ."^ 


1 “V n (3.1-3.31) 
'St ^ 
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From Eq. (3, 1.3. 6) 





^ 3.1.3.32) 


An additional useful relationship derived from Eq. (3.1.3.26) 
defines the difference between the inner ring and shaft 
displacement at their common surface, as a result of their 
rotational speed, in terms of the four integration constants 
of Eq. (3.1.3.28-31). 


O \o ^ S 

Under press fit conditions <S -O . 


(3.1.3.33) 


We will now examine the application of Eq. (3.1.3.28-31) 
in determining the ring behavior as a function of rotational 
speed. 


The following conditions might be encountered. 

1. A tight fit remains tight 

2. A tight fit loosens 

3. A loose fit remains loose 

4. A loose fit tightens 

For all four conditions, ® and rx* « O. ^ i.e. it is 

assumed that no change in pressure occurs at the internal 
surface of the shaft or the external surface of the inner 
ring resulting from ring rotation. Also Jlx and.H.« are 
identical. All four of the integration constants Ki* < . 

Kis and Kxs are unknowns. Also, either the change in pressure 
■pXL rotational speed, iX embodied inM^ and % are 

additional unknowns. 

Now the formulations of Eq. (3.1.3.28-33) are presented 
which are required to solve each of the four conditions. Since 
the equations refer to section SI, the variable t^jt^^will be 
written as U. 
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3. 1.3. 6.1 Situation 1 and 2. Initially Tight Fit 

o * 

In addition to the integration constants, Pit is the 
unknown, where is the change in inner rijig internal pressure 
resulting from the rotational speed. If PxC is less than 
from Eq. (3.1.3.23), situation 1 is realized. The operating 
fit pressure is the difference between and 

fit remains tight. The resulting radial displacement at 
the inner ring roller path is given by: 




'XL -- 






^%T. 

T2. 


( 3.1.3.34) 


and the total change in roller path radius at section SI is 
given by: 

= VJ V vJt A Ujl. (3.1.3.35 

is the change in rolling path radius which occurs as a 
result of a temperature change from 68®F, 

If is greater than Pxi t he solution to the problem 

requires several additional steps. Eq. (3.1.3.28), must be 
re-solved to determine the rotational speed, XL -ru at which 
the change in fit pressure is just equal to the initial fit 
pressure and . This rotational speed is termed the 

tight fit speed limit and is subscripted with the letter (TL). 
U.sing the integration constants thus determined and 
Eq. (3.1.3.26)j Ufi. ®ay be determined. 

Note, the subscripts in parenthesis should be interpreted 
as follows: The first subscript refers to the state of the 

initial fit (T-tight) or(L-Loose), the second subscript refers 
to the rotational speeds (T-Total) or (L-Lirait). For 
instance, refers to the change in inner ring rolling path 

radius due to the rotational speed at which the initially 
tight shaft fit is lost. 
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After determining Utu Eqs. (3.1.3.28-31) are resolved twice 
for the four integration constants, with Pxi set to zero. 

First, for the full rotational speeds, yielding Ut.T and then 
for the tight fit speed limit, yielding . This is 

accomplished by change in the values of Wt. and Hi) through a 
change in the value of in Eq. (3.1.3.25). Utilizing 

super-position, the change in rolling path radius resulting 
from rotation is given by Eq. (3.1.3.36). 



(3.1.3.26) 


Eq. (3.1.3.26) is presented graphically in Fig. 12. Then 
is calculated from Eq. (3.1.3,25). 


3. 1.3.6. 2 Situations 3 and 4 Initially Loose Fit 


The integration constants Kix , Kji , , and are 

unknowns. Px;"** . , P* I** and are all set to zero, and 

only Eqs. (3,1.3,28-31) are solved. Using the integration 
constants thus obtained, Eq. (3.1.3.33) is solved for , If 
5 *** <> 1 P«OT I where Pmot is negative, the loose fit remains 
loose at operating conditions. The constants Kix and Kax are 
used to calculate Uxt. . Using Eq. (3.1.3,25), vJtot' is 
calculated. If the shaft has expanded further 

under the effects of rotation than the inner ring, and the 
initially loose fit becomes tight. This requires additional 
solutions to Eqs. (3.1.3.28-33). 

C is set equal to Pmot and the equations are solved 
for the integration constants plus the speed at which the 
initially loose fit becomes tight Using the integration 

constants K|* and along with Eq. (3.1.3.37), Ull is 

determined. Eq. ( 3.1.3.28-33) are again resolved twice at XJ-tu 
andil-uTfor Pxi*' and the integration constants, after setting S’^ 
to zero. U-rt. and can be obtained by from Eq, (3.1.3.34) 

using the two sets of integration constants. The operating 
fit pressure is then given by; 

1^x1 - TT - rv- (3.1,3.37) 



(3.1.3.38) 
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and Eq. (3.1.3.35) is used to determine Utot . Where Ut©t 
represents the total change in the radius to the rolling path 
of section Si as the bearing is taken off the shelf at 68°F, 
mounted on the shaft subjected to load and rotated up to 
operating speed. 

3. 1.3. 7 Bearing Geometry Changes 

The preceeding sections dealt primarily with section SI. 
However, the equations presented are valid for all four of 
the assembly sections. For each section a value is determined 
for UfOT , cold fit pressure, operating fit pressure and the 
speed at which an initially tight fit loosens. The changes 
in bearing geometry are calculated by compiling the Utot terms 
from all four sections. 

As noted earlier, Eq. (3. 1,3. 9) can be used to determine 
this change in rolling element diameter such that: 

Urt 

The change in bearing diametral clearance is given by 


(3.1.3.39) 


(3.1.3.40) 


where TaT is omitted from the subscript. The change in contact 
angle at the outer raceway is given by: 




(3.1.3.41) 


At the inner raceway: 
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3.1.4 Bearing Fatigue Life 


Bearing raceway fatigue life, at a given laminutu as 

determined by Lundber g-P almgren , ref. 21 is expressed by 
Eq. (3. 1.4.1). 

4 

^ ' (3. 1.4.1) 


■ lO mu. 


QCr 


1_ Qe«iK 


Q CfAK. is the raceway laminum dynamic capacity, defined 
as the load for which the laminum will have a 90 percent 
assurance of surviving 1 million revolutions. m refers to raceway, 
k refers to laminum, from ref. 21 






(t ± cos 


1.074074 J 


o ,'7Tn?9> 

(3. 1.4. 2) 


.0.1Z22Z 

(y / [_S*( I i 


~I0.2S 


The upper sign is used for the outer race, the lower sign refers to 
the inner race. 


is raceway laminum equivalent load. 


Q 




.i- I 


Ve 


(3. 1.4. 3) 


is the individual roller contact load. On the k-th 
laminum and 4. = 4.0 or 4 = 4.5 depending respectively upon 

whether the applied load rotates or is stationary with respect 
to the raceway in question. 


Substitution of Eq. (3. 1.4. 3) into 


1 


L 




ft\K 


I 




cs?„, 










Eq. (3. 1.4.1) 

I 


yields 

(3. 1.4. 4) 


Since a unique lubricant film thickness can be determined 
for each roller/raceway contact, a unique lubricant life 
multiplier a^ is inserted into Eq. (3,1.4i4) to yield 
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L 


>omK 


^ ^ mxT y 

. 5 “. ^ Qcr.^ ) 


Q 


4/6 


V' \O^CCj 


(^30.-) 




(3. 1.4. 5) 


The Sg factor is a function of the EHD contact lubricant 
film thickness. Tallian, ref. 22 and Skuaka, ref. 23 , have 

defined the lubrication-life effect in terms of the lubricant 
film thickness/surface roughness ratio. 


A - ^ 


(3. 1.4.6) 


h is the lubricant film thickness and O" is the composite 
raceway roller surface roughness. 


J ^ 


r<\ 


1 


(3. 1.4. 7) 


The fatigue life data presented herein use the mean between 
the Tallian-Skurka data as shown by Fig. 13. 

The Lj^q life of a raceway is given by 

L 


lo 




RilJ 




I 


(3. 1.4.8) 


32 is a material factor. Background into the use of the 
material factor is presented in ref. 24 . For ease of comparison 

the life data presented herein, represent an a 2 factor of 1.0. 


The bearing fatigue life is given by 

r 


L 


\o 


rr\: X 


' ® A> 

( L Ibtfcfc7 

j _ — ^ 


(3. 1.4.9) 


Ni 


RPm 


N is the absolute value of the difference between the inner and 
outer ring speeds in rpm. 
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3.2 Analytical Procedure 

Presently, the greatest value of computer program AE71Y001 
is its ability to compare the relative performance of several 
bearing design or lubrication schemes as the bearing operates 
under identical load and speed conditions. The program was used 
in this manner in the bearing optimization phase of this effort. 

Use of S D§ [F Computer Program AE71Y001 requires prior de- 
finition of the following conditions. 

. Bearing design geometry and shaft housing dimensions 

. Applied load and/or outer ring angular misalignment 

. Shaft and/or housing rotational speed 

. Lubricant type, flow rate, input temperature, density, 
and tabulated viscosity temperature and pressure data 

. Ambient thermal conditions. 

3.2.1 Bearing Optimization Study Objectives 

(1) Based upon successful performance reported in the 
literature the maximum roller end-flange Hertz stress 
must not exceed 3.11 x 10*^ N/M , (45,000 psi) 

(2) For the given bearing section height and width, fatigue 
life is maximized to the extent that objective (1) 

all ows . 

(3) Bearing heat generation rates must not exceed the bear- 
ing and lubrication system's ability to dissipate the 
heat while maintaining a maximum lubricant temperature 
of 422°K (300°F) and/or a maximum component temperature 
of 450°K (350°F). 

The bearing design was optimized for the maximum NASA sup- 
plied load conditions as presented in Enclosure 1. 

. Axial load 53378 N (12000 lb) 

. Radial load 26689 N (6000 lb) 

Shaft speed 12500 rpm 5 Y 
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. Lubricant input temperature 340°K (195°F) 

. Maximum lubricant flow rate 0.017 (4.5 gal/min) 

min 

3.2.2 Roller End-Flange Contact (REFC) Performance Limitations 

Computer program AE71Y001 generates the following data which 
describes the REFC performance for each roller: 

. Contact load 

. Maximum Hertz stress 

. EHD film thickness 

. Coefficient of friction 

. Heat generation rate 

. Rolling and sliding speeds 

. Lubricant, flange and roller temperature. 

It remained to determine operational limits of the above 
variables which permit acceptable performance. 

Data on sliding speed and smearing load was presented by 
Cocks et. al. Ref. (25). This data, when extrapolated to 
the conditions of lower stress and higher temperature and sliding 
speeds, applicable in the phase of study, indicated that smearing 
would not occur as long as lubricant supply was adequate. 

More recently, Conners et. al. ref. (5) presented data on 
25 mm bore tapered roller bearings which were operated at REFC 
Hertz stresses to 2.69 x lo" N/m^ (39,000 psi) and at a flange 
linear speeds to 148 m/sec (3136 oft/min). Data from ref. (5) 
indicated that satisfactory performance could be obtained at 
the values of Hertz stress and sliding speed, required as long 
as the lubricant is readily available to the contacts. 

Lemanski, ref. (6), has provided additional data with which 
roller end flange contact performance could be evaluated. The data 
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is for a 95.3 mm (3.75 inch) bore bearing operating at 16000 rpm, 
resulting in a flange linear speed of 101.6m/sec (20,000 ft/min) 
versus 107 m/sec (21,000 ft/min) for the bearing proposed. Hertz 
stress values ranged from 1.93 x 10^ to 3.65 x lo" N/m2 
(28,000 to 53,000 psi) which corresponds very closely to the 
stress levels encountered in this study. A MIL-L-7808G 
lubricant was used in the ref., (6) experimental work as compared 
to the more viscous MIL-L-23699 used in this analysis. 

Lemanski's experimental results at a flange linear speed 
of 101.6m/sec indicate successful REFC performance at the 
following conditions. 

9 9 9 

Hertz stresses from 2.07 x 10 to 3.58 x 10 N/m*^ 

(30,000 psi to 52,000 psi). 

. Lubricant flow rates of 98 Kg/Hour (215 Ib/hour) 

introduced through the race to the roller large end 
flange and 47 Kg/hr (105 Ib/hr) introduced by jet at 
the roller small end side. 

The data from all three sources were used to establish a 
basis to Judge the roller end flange performance at the 
various contact angles examined in the optimization study. 

In particular, a maximum Hertz stress level of 3.1 3 x 10^ N/m2 
(45,000 psi) was used as the limiting criteria. Sliding and 
flange linear speeds encountered in this study have been 
successfully exceeded in the referenced tests. 

3.2.3 Phase 1 - Preliminary Design Considerations 

At the outset of the design study it was anticipated that 
outer ring contact angle would have two prime effects on bearing 
performance which would have to be considered in selecting an 
optimum bearing design. 

1) With increasing contact angle bearing raceway fatigue life 
was expected to increase. 

2) The operating conditions at the roller end flange contact 
(REFC) were expected to become more severe as contact 
angle increased. 
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Fatigue life was expected to increase with contact angle 
due to the two to one ratio of the applied axial to radial load 
and the reduced roller race normal loads which would result at 
high-contact angles. The operating conditions at the REFC were 
expected to become more severe. As contact angle increases the 
bearing apex length decreases; the REFC principal radii decrease 
and as a result the Hertz stress increases. Outer ring contact 
angle the number of rollers and roller dimensions had to be 
established which would fit within the limits of bearing radial 
section height and width and in addition, would provide maximum 
load carrying capacity. An analysis was developed for this 
purpose and is presented in Appendix A. 

Using the methods of appendix A a unique set of roller 
dimensions, large end diameter, length included angle and 
number of rollers were determined for various contact angles 
ranging from 12° to 28°. This data is presented in Table 1. 

3.2.4 Outer Rina Contact Angle Determination 

Using computer program AE71Y001 and the bearing designs 
as determined by the methods of Appendix A, bearing performance 
analyses were made for outer raceway contact angles of 12, 16, 18, 

20, 22, 24 and 28 degrees. Fig. (14-17) present the relevant 
data with which the bearing contact angle was selected. This 
data was generated using a roller sphere end radius of 90 percent 
of the apex length for all designs and a constant lubricant 
temperature of 364°K (195°F). 

From Fig. 14, bearing fatigue life increases with increasing 
contact angle due to the decrease in roller raceway normal loads, 
apparent from Fig, 15. 

Fig. 16 presents a relatively complete picture of the 
roller end contact performance. REFC Hertz stress increases 
rapidly with increasing contact angle. Film thickness remains 
fairly constant. Fig, 17 presents the bearing frictional 
heat generation rates which increase with increasing contact 
angle . 

Utilizing the data of Figs. (14-17) nominal outer raceway 
contact angle of 24° was selected as the best compromise between 
maximizing fatigue life and maintaining an REFC Hertz stress level 

below 3.11 X lo" N/m^ (45000 psi), 
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3.2.5 Cage Design 

As a part of this design study, consideration of an optimum 
lubrication scheme played a prime role. A large quantity of 
lubricant within the bearing cavity serves only to increase 
frictional heat generation and to cause increased slip in the 
roller raceway contacts by increasing the drag force on the 
rollers . 

Efficient lubrication necessitates getting the lubricant 
to the bearing contacts which must be lubricated and then 
removing the lubricant as quickly as possible. The "Z" cage, 

Fig, 0, was selected because it does not impede efficient 
lubrication. 

The "Z" cage is guided by the inner ring land at the roller 
large end side of the bearing and by the outer race at the small 
end side. This configuration permits relatively free entrance 
of Jetted lubricant at the small end between inner ring land 
and cage. Since centrifugal force will drive the lubricant to 
the outer raceway once it is in the bearing, it is certainly 
beneficial to aim the jet at the inner raceway. Centrifugal 
force and the outer raceway taper tend to drive the lubricant 
to the large end side of the outer race. As a result, 
undesirable lubricant churning may be minimized if the lubricant 
is free to exit from the bearing at that point. The "Z" 
cage allows this to occur since the cage is inner ring land riding 
at the large end. 

3.2.6 Bearing Manufacturing Considerations 

The basic roller configuration as presented in Table I 
for the 24° outer raceway contact angle was reviewed by bearing 
design specialists who finalized the design, making certain that 
the bearing could be manufactured and that it would conform to 
accepted design principles. It was stipulated that the cage be 
of the "Z" configuration. The following design changes resulted 
from this review. 
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1) The roller length was reduced from 33.934 to 27.742 mm 
(1.3360 to 1.0922 inches) 

2) The roller large end diameter was increased from 23,564 
to 23.609 mm (0.9277 to 0.9295 inches) 

33 The roller included angle was increased from 5° 36’ to 
5° 40' 

4) The number of rollers was increased from 19 to 20. 

The "Z" cage was primarily responsible for the most important 
of these design changes, the reduction of roller length. For 
this reason, an inner ring land riding cage was examined. 

As noted previously, the Z cage was chosen for its ability 
to permit easy entrance of a jet stream between the cage rail 
I.D. and the inner ring land O.D. at the small end side of the 
bearing. An inner ring land riding cage severely limits the 
jet target and perhaps could cause a lubricant starved condition 
at the roller small end side of the bearing. Two alternatives 
exist to circumvent this situation; 

1) Introduce lubricant through the race to the undercut 
at the roller small end side of the inner race. 

2) To scavenge the lubricant only at the large end side. 

The latter may result in i nef f i cient 1 ubr i cat i on, and 
thus was not considered further, 

A design review considering through race lubrication to 
the small end undercut, revealed that undercuts at both the 
large and small ends of the inner race would have to be 
substantially over sized in order to guarantee proper placement 
of the lubrication holes. This situation was caused by the 
large tolerance stack up alloted to ring machining dimensions. 

The decrease in inner raceway effective length resulting from 
the large undercuts approached the reduction in roller length 
required by the "Z" cage. The "Z" cage was thus retained in the 
final design. 
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3.2„7 Roller End Flange Contact Lubrication and Geometry 

The roller end flange geometry is of substantial importance 
in the high speed tapered roller bearing. Since that contact 
has been shown to be most prone to failure. 

The data and conclusions of ref. (5, 6 and 25) all support 
the need for a copious supply of lubricant at the flange roller 
end contact. For this reason through race lubrication to the 
undercut at the roller large end side of the inner race is 
specified for the proposed optimum design. The through race 
flow hole di ameter and spacing is in accordance with the methods 
of ref. (6). The amount of lubricant to be supplied will be 
discussed subsequently. 

Also of prime importance in maximizing the chances of 
survival of the REFC is the roller end geometry. Historically 
a roller end sphere radius between 95 and 100 percent of 

the bearing apex length has been used. However, Lemanski 
encountered problems at high speed operation using an Rgpjj 
equal to 97 percent of the Apex length and found optimum results 
with Fgpfj equal to 80 percent of Apex. Conners used an Rgpjj 
equal to 75 percent of Apex in his tests. 

Using Archard's EHD film thickness equation, Lemanski 
plotted REFC film thickness versus K where: 

V _ 3.2.1 

tv- /ypELX 


He found maximum film thickness between K = 0.75 and 0.80. 


Archard developed a term 4^ in his film thickness equation 
to express the effect of side leakage in an elliptical contact. 

4 = t 4 - 3.2.2 

L 3R^J 

where y j^s the direction of translation. For a typical line 
contact is infinite and 4^ = 1 . In his theoretical derivation 

Archard applied an exponent on of 0.74 in his film thickness 
equation. However, in his experimental work he found that an 
exponent of 0.93 yielded a better fit of theoretical to 
experimental data. As indicated by Eq, (3.1.1.26) the exponent 
0.93 has been used in this anal^^is. 
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Lemanski’s determination of a maximum film thickness at 
K b^ween 0.75 and 0.80 can be explained by expressing 
and"^^, Eq. (3.1.1.28 and 3,1,1.29) in terms of K°Apex and 
then solving the expression for K, h is the lubricant 

film thickness. Using the exponent of 0.74, K is found to be 
0.754. With the ^ exponent of 0.93, K is found to be 0.66, 

Lemanski's findings suggest experimental justification for 
use of a K value less than the historical value of 0,95 
to 1.00. Archard's film thickness equation also suggests 
justification. 

Computer Program AE71Y001 was used to obtain data on the 
roller end flange contact performance characteristics as a 
function of K for the proposed optimum bearing design. These 
data are presented in Fig. 18. Based on this data and the 
experimental results of ref. (5) and (6) the roller end sphere 
radius equal to 82,5 + 2,5 percent of the Apex length was 
chosen for the proposed bearing design. 

3.2.8 Optimum Tapered Roller Bearing Thermal Performances 

Of prime interest in this effort was the development of a 
model to predict the thermal performance of the tapered roller 
bearing including the effects of the lubrication input and 
scavenge system. 

The bearing-shaft-housing assembly dimensions used in the 
thermal portion of the analysis are shown in Fig. 19 , The 
thermal boundary and lubrication conditions studied are shown 
in Fig. 20, A schematic of the lubrication distribution 
schemes studied is presented in Fig. 21. 

The thermal analysis section 3.1.2 was used to evaluate the 
various lubricant flow rates to the bearing' by jet to the roller 
small end; and by thicugh race lubrication to the large end 
undercut. In generating these data the assumption was made that 
the amount of lubricant within the bearing cavity could be 
represented as a function of total lubricant flow rate 
according to Fq. (3.2.3) 

■P = Z-^.O '<lo^ M 3*2.3 
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^Where P is the decimal fraction of lubricant in the cavity 
and M is the total lubricant input rate to the bearing. Normally 
P and M are separate inputs to the computer program. This 
assumption is believed to be conservative, i.e. the heat 
generation rate predictions which result from the assumption 
will be high rather than low. Thus, if the assumption is not 
totally accurate, less rather than more lubricant should be 
required to maintain desirable system temperatures. 

The proposed bearing lubrication scheme consists of 
lubricant input to the undercut at the large end side of the 
inner race and by jets to the small end side of the bearing 
between the cage rail I,D, and the inner ring land O.D, The 
lubricant is scavenged at both sides of the bearing. With 
this lubrication scheme, an optimum amount and ratio of through 
race to jet lubrication was sought. Figs. 22-29 present the heat 
generation rate and lubricant temperatures at the flange, inner 
raceway, outer raceway and the average lubricant temperature 
which results from the various flow rates and flow ratios. 

Figs. (21-25) indicate that a jet lubricant flow rate of 
approximately 340 kg/hr (750 Ib/hr) is required to maintain 
the majority of the lubricant temperatures between 394 and 
422°K (250 and 300®F), Only 66.2 percent of the lubricant 
jetted to the small end side is calculated to enter the 
bearing. 33.8 percent is rejected to the scavenge. The lubricant 
input through the race to the large end undercut effectively cools 
the flange. Fig. 19. However, it is not well distributed 
throughout the bearing. Fig. 22 indicates that a through race 
flow rate of 227 kg/hr (500 Ib/hr) is sufficient to maintain 
the temperature of the lubricant at the flange to a desirable 
394°K (250°F). The shape of the curves representing the 
temperature versus through race flow rate in Fig. 22-25 indicate 
that the cooling ability of the through race lubrication tends 
to become nominal above a flow rate of 227 Kg/hr. Referring to 
Fig. 21.2, the lubricant distribution model assumes that the 
major portion of the lubricant enters at the undercut and moves 
to the outer race node 3 and then exits from the bearing. Only 
nominal amounts of lubricant are assumed to be distributed to 
nodes, 4, 5 13 and 14. where frictional heat is being generated. 
Referring to Fig. 21.1 only 66,2 percent of the lubricant which 
is jetted to the bearing is assumed to enter. However, that 
which does enter is assumed to be better distributed over the 

inner race than when the lubricant enters through the race. 
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Fig, 23-24 indicate that despite less lubricant actually in 
the bearing, a better distribution of the lubricant results in 
lower lubricant temperatures at inner and outer raceways with 
jet as opposed to through race lubrication at flow rates greater 
than 272 kg/hr 600 (Ib/hr). 

The data presented in Figs. 26-29 was generated to determine 
the effectiveness of introducing the lubricant via. jets at the 
small end side of the bearing, versus introducing it directly to 
the large end flange. Difference in lubricant temperatures at 
the same total flow rate was used as the measure of effectiveness. 
Each curve has been generated by holding one of the lubricant 
input rates constant at 227 kg/hr (500 Ib/hr) while varying the 
other from 56,7 to 227 kg/hr (125 to 500 Ib/hr), in increments 
of 56,7 kg/hr. Along any line parallel to the ordinate the 
total lubricant flow rate is identical, ranging from 283.5 
to 453.6 kg/hr (625 to 1000 Ib/hr). The difference in temperature 
from curve to curve is thus a function only of the method used 
to introduce the lubricant. 

Figs. 27 and 28 show lubricant tepperatures at the small 
end of the inner and outer raceways as well as the center and 
indicate lower temperature at the small end side of the raceways 
resulting with increased jet lubrication. 

With the beneficial effects realized at the small side 
of the raceway, the relatively uniform and low temperatures 
throughout the bearing cavity, and the requirement for lubricant 
at the small end of the outer raceway for cage lubrication, 
lubricant flow rates of 227 kg/hr (500 Ib/hr) via. jet 
and through the race were selected for optimum lubrication. 

3,2.9 Geometry Change at Operating Conditions 

In the optimization study, the geometry change analysis 
was used to obtain changes in b earing geometry which would be 
manifested at operating conditions. 

A cold shaft fit of 0.114 mm (0.0045 in) tight and a 
line to line housing fit were assumed. With the recommended 
lubrication scheme and flow rate?, at the heavy load condition, 
two significant changes in geometry were noted. 
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1} Changes in raceway and roller diameters result in a 
decrease in bearing diametral clearance of 0.089 mm 
(0.0035 in). More appropriate to tapered roller 
bearings, this corresponds to a reduction in shaft 
end play of 0.1119 mm (0.0044 in). 

2) The outer raceway contact angle is expecte to increase 
by 2.2 minutes and the inner raceway contact angle 
is expected to decrease by 0.2 minutes. 

Sufficient bearing end play must be built into the bearing, 
shaft housing assembly such that at operating conditions the 
bearings do not become axially preloaded. 

The indicated changes in bearing contact angle would result 
in a pinching across the small end of the roller. This condition 
can be alleviated by the bearing design. The optimum bearing 
design thus has a specified outer raceway contact angle 2,4 minutes 
less than nominal. 

3.3 Analytical Results 

Since, in the bearing optimization s t udy, r esul t s of one 
phase were needed to proceed to the next phase a good portion 
of the analytical results have been presented in the Analytical 
Procedure Section. However, those results not previously noted 
are presented herein. 

3.3.1 Bearing Material, Lubricant Film Thickness and 

Bearing Fatigue Life 

As noted, in Section 3.1.4 on fatigue life analysis, a 
multiplicative factor is applied tc raceway fatigue life to 
include the effects of material and lubrication. Throughout, 
the material factor was set to one. However, as noted, the 
effect of lubrication as determined by the.^. ratio, varied with 
EHD film thicknesses. In the preliminary studies in which an 
optimum contact angle was being sought, the lubricant temperature 
was set to 364°K (195°F). The resultant ^nner and outer raceway 
film thicknesses were typically 127 x 10 mm (50 microinches) 
with a ratio of approximately 7, and the life multiplier was 
2.85. In the later studies in which the lubricant temperature 
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was approximately 394°K (250°F), the film thickness dropped to 
(68.6 X 10“^ m) 27 microinches yielding a ratio of 3.85. 

However, the life multiplier dropped only to 2.58. This data 
indicates that nearly full EHD films exist in the roller raceway 
contacts with the proposed design at the proposed lubricant 
flow rates. Also, a substantial increase in bearing life can 
not be expected with increasing lubricant film thickness beyond 
approximately 76.2 x 10~^ m (30 microinches). 

3.3.2 Contact Load and Cage Slip 

The effect of raceway normal loads on cage slip can be noted 
in Fig, 15. The high contact loads which occur at the 12° 
contact angle result in the cage speed deviating from epicyclic by 
less than 1 percent. As the raceway normal loads decrease with 
increasing contact angle, cage speed deviation from epicyclic is 
6,9 percent at 28°, Cage speed deviation from epicyclic (cage 
slip) is indicative of slip within the roller inner raceway 
contacts. This slip explains the increase in contact heat 
generation rates with increasing contact angle evident in Fig, 17. 

The decrease in churning heat generation with increasing 
contact angle simply reflects the decrease in roller frontal area. 

The effect of reduced contact loads on cage slip can also 
be seen in Fig, 30. Fig. 30 is a plot of percentage cage slip 
versus applied radial load for the proposed bearing. The data 
of Fig. 30 was generated in order that the performance of the 
optimum bearing design could be examined at the low and midpoint 
of the load spectrum with a lubricant flow rate of 454Kg/Hr 
(1000 Ib/Hr), half entering the bearing via. the small end jets 
and half through the inner ring to the large end undercut, the 
relatively light loading results in substantial skidding in the 
inner raceway-roller contacts. The high slip rates under the 
low load-high speed condition are indicative of too much lubricant 
within the bearing cavity. This finding indicates 

(1) That in order to have the bearing run at an acceptable 
slip rate, (less than 10 percent), the lubricant supply 
will have to be reduced with possible higher temperatures 
resulting. 
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(2) The relationship used to calculate roller drag Eq, , 
(3.1.1.34) is not totally accurate or the assumed 
relationship between lubricant flow rate and the amount 
of lubricant with the cavity Eq. (3.2.3) is not totally 
accurate. 

Most likely, both situations (1) and (2) are true. 

The bearing lives at the light and medium load conditions 
are presented in Fig. 31. These lives include the effect of 
reduced roller centrifugal force caused by the high slip rates. 
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4.0 DISCUSSION OF RESULTS 


Two prime results have emerged from this effort. The capabilities 
of computer program AE71Y001 have been enhanced and refined such 
that it is now a powerful tool for the kinematic and thermal 
performance analysis of high speed tapered roller bearings. 

The second prime result is the proposed optimum tapered 
roller bearing design, developed by means of the computer program. 

The complete bearing design is presented in Enclosure 2. The 
performance characteristics for the bearing with the proposed 
optimum design and lubrication scheme is presented in Table II. 

The plots in Fig. 14-18 which were used in determining the 
optimum bearing contact angles are discontinuous at an angle 
between 20° and 22°. Appendix A, which was used to determine the 
roller geometry at each contact angle, contains three criteria to 
limit roller size; tests T^, T 2 and T 3 . At a given contact angle, 
the roller fully utilizes the dimensions of only two of the 
tests. Roller sixe was limited at outer raceway contact angles 
by tests T^ and T 2 , below and including 20°. Tests Tj^ and Tg 
limited size at angles above and including 22°. This change 
in criteria appears to be the cause for the discontinuities. 

It became apparent from the study that both on the basis of 
a bearing's load carrying capability and lubrication requirements, 
a bearing design which would perform optimally at high load 
conditions, would not necessarily be the design chosen at lower 
load conditions, especially if the ratio between axial and 
radial loads were different at the high and low load conditions. 

The bearing design optimization was directed only at the high end 
of the load spectrum presented in Enclosure 1. It was reasoned 
that any bearing which would perform under the most severe 
load conditions would also perform satisfactorily at lower load 
conditions. The inverse is probably not true. 

Computer program AE71Y001 is capable of analyzing two 
identical bearings mounted on the same shaft, such that an 
axial load in either direction will be taken by one or the other 
bearing. This study, as required, considered only a single 
bearing. Use of the data in Fig, 30 and 31 was generated using 
applied radial load only, since the axial loads required to keep 
the bearing together under the radial load were unknown. In 
general, only a two bearing system is capable of accepting applied 
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loading in five distinct directions. Computer program AETiyoOI 
is capable of analyzing the two bearing assembly, but this 
capability was not used in this study. 

As noted in the introduction, the principle value of computer 
program AE71Y001 is its ability to compare various bearing 
designs and lubrication schemes as the bearing operates under 
identical loading and operating speed conditions. The program's 
ability to predict absolute magnitudes of bearing performance 
parameters particularly in the area of heat generation rates 
and system temperatures must be evaluated as experimental data 
becomes available. The lubrication scheme and lubricant flow 
rates recommended will be only as accurate as the assumptions and 
relationships used in determining heat generation rates and the 
lubricant distribution. However, as mentioned, the recommended 
flow rates are believed to be conservat i ve . Even though 
conservative, the recommended flow rates are well within the 
system limits as specified in Enclosure 1. 
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5.0 CONCLUDING REMARKS AND PROPOSED FUTURE EFFORTS 

Consistent with the objectives of this effort an analytical 
method has been developed to predict the thermal and kinetic 
performance of lubricated tapered roller bearings subject to 
combined loading. This analysis, D§ [F computer program AE71Y001, 
provides the user with a powerful tool to analyze new bearings 
designs (as shown by the bearing optimization study) and to 
trouble-shoot problem applications. 

5.1 Analytical and Experimental Data Correlation 

Effort should now be made to compare the program predictions 
with experimental data. Hopefully an effort in data correlation 
would result in an analytical expression relating bearing size, 
speed, and lubricant flow rates to 1 ubri cant -imposed resistance 
to roller motion. A technical break-through in this area is 
needed and if found could be applied to other bearing types. 

A second but as important result of a data match-up program 
would be knowledge gained of the strong and weak areas of the 
mathematical model. Weak areas could be subjected to additional 
analysis and thus strengthened. 

5.2. Study of Concentrated Contacts 

Another area for future effort is the study of low stress, 
high sliding rate-cont acts such as occur in the roller end-flange 
contact (REFC). Ideally the objective of this study would be 
to analytically and experimentally deduce a performance number 
such as; . 

-Njr-- Sr(lA, ,5.1) 

h = lubricant film thickness 
H = maximum Hertz stress 

V = sliding rate 

u = lubricant entrainment velocity* 

• 

M = lubricant supply rate 

= Archard's ellipticity. factor 
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T = temperature 

Based on experimental results limits of could be established 
which would denote excellent, adequate and poor performance. Such 
a development is essential to the rational design of low stress, 
high sliding rate-concentrated contacts. 
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TABLE I - PRELIMINARY BEARING DESIGNS 


Contact 

Angle 

(Dea) 

^ Roller 
Included 
Angle (Deg) 

Roller 

Length 

(MM) 

Roller 

Length 

(in) 

Roller 
Large End 
Diameter (MM) 

Roller 
Large End 
Diameter (in) 

Apex 

Length 

(MM) 

Apex 

Length 

(in) 

Number 

of 

Rol lers 

12 

1.70 

32.994 

1.2990 

27.960 

1 . 1008 

471.22 

18.552 

15 

18 

2.33 

33.538 

1.3204 

25.798 

1.0157 

317.29 

12.492 

17 

20 

2.51 

33.660 

1.3252 

25.121 

.9890 

286.79 

11.291 

18 

22 

2.67 

33.667 

1.3255 

24.397 

.9605 

261.87 

10.310 

18 

24 

2.80 

33.934 

1.3360 

23.563 

.9277 

241.15 

9.494 

19 

28 

3.04 

33.688 

1.3263 

22.164 

.8726 

208.96 

8.227 

21 


z 
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TABLE 2 


PERFORMANCE PARAMETERS AT MAXIMUM LOAD CONDITIONS FOR 
THE PROPOSED OPTIMUM TAPERED ROLLER BEARING 
SEE ENCLOSURE 2 FOR DESIGN DETAILS 



Max 

N 

. Contact 
Load 

LBF 

Maximum Hertz 
Stress 

10^N/m2 lO^Dsi 

Minimum Lub. 
Film 

Microns Microin. 

Sliding Heat 
Generation 
Watts BTU/Hr 

Li?s 

Hours 

Surface Temp. 
- °K °F 

Outer Ring 

9525.8 

2141.5 

1 . 1059 

1.604 

.7004 

27.576 

95.49 

325.91 

1308. 

409.62 

277.67 

Inner Ring 

7088.6 

1593.6 

1.0397 

1.508 

.6865 

27.030 

140.51 

479.55 

1547. 

384.80 

233.00 

Flange 

1731.6 

389.3 

.3013 

.437 

.9773 

38.477 

92.51 

315.74 

- 

375.21 

215.73 

..j 

Max. 

Cage Load 

Rotational 

Gyroscopic Moment 

Centrifugal Force Orbital 

Skewi ng 

Moment 

a> 

N 

LBF 

Sneed 

RPM 

N-M 

IN-LB 

N 

LBF 

Sneed RPM N-M 

IN-LB 

Roller 

85.02 

19,113 

47450 

8.342 

73.840 

2536.9 

570.33 

5146.9 

.1422 

1.2586 



Linear Deflection 

Linear Snrina Rate 



LIO 

Cage 


X-Axis Z-Axis 

10"^ lO'^IN 10~^M lO'^lN 

DFX/DX, , DFZ/DZ 

10"^N/M 10°LB/IN 10°N/M lO^LB/lN 

Torque 

N-M IN-LB 

Li fe 
Hour s 

Slip 
2 

Bearing 

.047 1.882 .0081 .3198 

.6244 5.527 2.105 18.628 

8.161 

72.235 

765.2 

5.774 


Heat Generation Rate 
Watts BTU/Hr 

10688. 36478. 


Shaft Speed 12500 rpm 

Total lubricant flow rate 454Kg/hr ( 1000 Ib/hr) 

Lubricant inlet temperature 364 °K (195°F) 

Lubricant type MIL-L-23699 

Roller sphere end radius is 80 percent of the apex length 
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Fa 



FIG. 1 Bearing operating in YZ - Plane 
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0.62 X 13.75 
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FIG. 3> ROLLER RACEWAY STRESS CONCENTRATIONS 



(b) PARTIALLY CROWNED ROLLER STRESS CONCENTRATIONS 


-K- EFFECTIVE LENGTH OF ROLLER - RACEWAY CONTACT 


I 
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R 


CROWN 



^LAT W 


W j 


/ 



LAMINATED ROLLER STRUCTURE SHOWING 
CROWN DROP VERSUS LENGTH 


F,6 4 TAPERED ROLLER GEOMETRY 
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FIG. s Components OF Roller Deflection Due To 
. Radial Load, Misalignment AMD CROWNING 
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ABSOLUTE VISCOSITY- ME WTO N SECOND/METER 


' • ' • ' ; . ALT3P010 

FIG. 6 Mil-L-23699 Lubricant Viscosity VS 
Pressure AND Temperature ' : 



F.4186A 
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FIG. 7 X2 PLANE INNER RING 

DISPLACEMENTS 
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FIG. B CAGE TY 


ROLLED RIDING 


INNER RING LAND RIDING 


OUTEA RING LAND RIDING 

INNER AND OUTER 

RING 

(ORR) 

RIDING ( Z. 
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^ 63 64 - 65 66 67 68 

Circles Represent Temperature Nodes in the 
Assembly Components 

SauAREs Represent Temperature Nodes inThe 
Lubricant And Fluid Coolant 
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FIG. lO.I LUBRICANT DISTRIBUTION MODEL 
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Designates the total mass flow rate from node 21 


Designates the rate at which lubricant flows fron node 21 
to node J 

is the scoop efficience; that fraction of the lubricant Jetted to the 
scoop from node 21 which actually enters the bearing, must be 
input to the analysis 

Is the flow rate from nodes and which enters the bearing 
vis the scoop flow paths. 

Is a function of bearing cage 

are auxiliary mass flow rate variables 

Is a function of shaft speed, the number of radial under race 
flow paths and flow path dimensions Eq. (3. 1.2. 2) 
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10.5 FuONaJ ScWEHE. FV 2 . 0 WI TvA-E. 
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Designates the total mass flow rate through node 11 

Designates the rate at which lubricant flows from node 11 to node J 
Is the jet efficiency that fraction of the lubricant jetted to 
the bearing. can be inpute or calculated as a function of jet 
diameter and bearing design 

js the flow rate from nodes 11 which enters the bearing from the jet 
are prescribed func ions of and Eq. (3.1.2.38.1 - 3.1.2.30.3) 
and are auxiliary mass flow rate variables 
is a function of bearing cage speed and Jet velocity Eq 3.1.2.26 
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FIG. ILZ BEARING ASSEMBLY 

WIDTHS, AND 
ROLLER PATH RADII 
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LIFE ADJUSTMENT FACTOR 
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FIG. 13 Life adjustment factor for Lubrication 0. 


Film thickness/surface Finish ratio A 



A 
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FIG. 15 ROLLER - RING NORMAL LOADS, ROLLER 

CENTRIFUGAL FORCE AND CAGE SLIP PERCENT 

VS 

OUTER RING CONTACT AfIGLE 


ROLLER END SPHERE RADIUS 
80% OF APEX LENGTH 

SEE TABLE- i FOR GEOMETRY DATA 



OUTER RING CONTACT ANGLE - DEGREES 
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PLANCE LOAP ~ HEW TONS (lBF) 
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FIG. FLANGE LOAD, HERTZ STRESS AND 
LUBRICANT FILM THICKNESS 

VS 


CONTACT ANGLE 

AXIAL LOAD* 53378 N (12000 LB) 
RADIAL L0AD = 26689N (6000 LB) 
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FRICTIONAL HEAT GENERATION RATES* WATTS (lO BTU/HR) 
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FIG. 1 7 FRICTIONAL HEAT GENERATION RATES 

VS 

CONTACT ANGLE 

SEE TABLE 1 FOR GEOMETRY DATA 


kso) 


14000-j 


130004 




12000-4 


IIOOOH 


100004 


k40) 




k^) 




RADIAL LOAD 26689 N feooo UB.) 
SHAFTSPEED 12500 RPM 



12 16 20 
CONTACT ANGLE 


24 2 6 

DEGREES 
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FIG. IS FLANGE HERTS STRESS; LUBRICANT FILM 
THICKNESS and HEAT GENERATION RATE 

VS 

ROLLER END SPHERE RADIUS 

CONTASr ANGLE = Z4-® 

AXIAL LOAD = 53378 N («Eooo Lb) 
RADIAL LOAD = 26669 N (sooo lb) 
SHAFT SPEED 12 500 RPM 



PERCENTAGE OF APEX LENGTH 
FOR ROLLER ENDSPHERE RADIUS 
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□ 

^ DIAMETERS, INCHES , WIDTHS OR EFFECTIVE LENGTHS, INCHES 


(0 

H 

A 

HOUSING O.D. 

m 

273.50 

Inch 

10-.750 

G 

HOUSING 

203.20 

Inch 

8.000 


B 

BEARING O.D. 

206.38 

8.125 

H 

OUTER RING 

38.10. 

1.500 

— 

C 

HOUSING SHOULDER DIAMETER 

190.50 

7,500 

I 

INNER RING ’ 

1*7.63 

. 1.875 

ni 

D 

SHAFT SHOULDER DIAMETER 

1W.05 

• 5.750 

J 

SHAFT 

165.10 

6.500 

(D 

E 

BEARING BORE DIAMETER 

120.65 

>1.750 

K 

OUTER RING-INNER RING 

9.53 

0.375 

z 

F 

SHAFT I.D. 

88,90 

3.500 


(CUP-CONE) OFFSET 



0 


FIG. 19 BEARING- 

•SHAFT 

AND HOUSING 

ASSEMBLY DIMENSIONS 

• 
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Scavenge Scavenge 



FIG. 20 BEARING BOUNDARY CONDITIONS 
AND LUBRICATION SCHEME 


Lubrication Data 

Lubricant Type MIL-L-23699 

Input Via Jets 

Number of Jets 2 

Jet, Diameter 2,54mm (0.10 in.) 

Flow Rate 227 KG/Hr (500 Ib/Hr) 

Jet Efficiency 0.662 

Inlet Temperature 264®K (195°F) 

Input Through Race 

. Number of Flow Holes 30 
Flow Hole Dia, 1.02ram 0.4°in. 

Flow Rate 227 Kg/Hr (500 Ib/Hr) 
Inlet Temperature 364°K (195®F) 

Lubricant Exit Vi'a 
.Scavenges at Both Sides 
Of the Bearing 
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^ Boundary Conditions 
Air ■ 

Temperature 293®K (75®F) 
Velocity - Quiescent 

Housing 

Material Steel 
End Temperatures 

Side 1 372®K (210®F) 
Side 2 364®K (195®F) 

Shaft 

Material Steel 
End Temperatures 

Side I 372®K (210®F) 
Side 2 364®K (195®F) 




Fig. 21 Lubrication Schemes Examined In The 
Bearing Optimizat|.on Study 



21.1 JET, At The Roller Small End To The 
Gap Between The Cage and Inner Ring 
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FRICTIONAL HEAT CENEPATION RATE - WATTS (bTU/HR X 10^) 
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FIG. 22 FLANGE LUBRICANT TEMPERATURE AND 
FRICTIONAL HEAT GENERATION RATE 

VS 

LUBRICANT FLOW RATE 
AXIAL LOAD 53378 N (izooolbf) 
RADIAL LOAD 26689 N (6000 LBf) 
SHAFT5PEED 12500 R PM 

SEE ENCL.2 FOR DESIGN DETAILS 
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( 200 ) ( 500 ) ( 750 ) ( 1000 ) 

1 * I ‘ 1 * ' , * 

(00 200 300 400 

LUBRICANT FLOW RATES - KG/HR (lb / hr) 
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FRICTIONAL HEAT GENERATION RATE -WATTS (btu/hr x io*j 
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PIG. 23 INNER RACEWAY LUBRICANT TEMPERATURE 
AND frictional HEAT GENERATION RATE 

VS 

LUBRICANT FLOW RATE 

AXIAL LOAD 53378 N (\ 2.000 lbp) 
RADIAL LOAD 26689 N (6000 lbf) 
SHAFT SPEED 12500 RPM 

^EE ENCLOSURE 2 POP DESIGN DETAILS 
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FIG. 24- OUTER RACEWAY LUBRICANT TEMPERATURE 
AND frictional HEAT GENERATION RATE 

LUBRICANT FLOW RATE 

AXIAL LOAD 53376 N (laoooLBF^ 
RADIAL LOAD £6669 N ('eooo Lbp) 

SHAFT SPEED l£500 RPM 



LUBRICANT FLOW RATES - KG/ HR (lbm/ hr) 


o 

o 

K 

< 

j ■ 105 

RESEARCH LABORATORY SICF* INDUSTRIES. INC. 



AL73P010 


FIG. 25 LUBRICANT MEAN TEMPERATURE 
AND FRICTIONAL HEAT GENERATION 
rate due to LUBRICANT CHURNING 

VS 

LUBRICANT MASS FLOW RATE 


AXIAL LOAD 5337e N (izodolb') 
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PIG. 26 FLANGE LUBRICANT TEMPERATURE AND 
FRICTIONAL HEAT GENERATION RATE 

VS 

LUBRICANT FLOW RATE 
SEE ENCLOSURE 2 FOR DESIGN DETAILS 

CONSTANT FLOW RATE= EZ6.79 KG/HR (sooLBM/ hr) 

AXIAL LOAD 5337 8 N (izooo LBf) 

^ RADIAL LOAD 266B9 N (sooo LBF) 

"o SHAFT SPEED laSOO RPM 



50 100 150 200 ZSO 


VARIABLE LUBRICANT FLOW RATE -KG/HR (lB / HR) 
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FIG. 27 INNER RACEWAY LUBRICANT TEMPERATURES 
AT CENTER AND SMALL END SIDE OF RACE 

and frictional heat generation rate 

LUBRICANT FLOW RATE 
CONSTANT FLOW RATES=226.79 KG/HR (sool&m/ha) 

AXIAL load 53378 N (l^ooo LBf) 
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FIG. 2a OUTER RACEWAY LUBRICANT TEMPERATURE 
AT CENTER AND SMALL END SIDE OFRACE 
AND FRICTIONAL HEAT GENERATION RATE 

VS 

LUBRICANT FLOW RAtE 

CONSTANT FLOW RATE = 226-79 KG/ HR 

CSOO LBM/HR) 

AXIAL LOAD 53378 N (izooo L0 f) 

RADIAL LOAD 26689 N (sooo lbf) 
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Fig.2® lubricant bulk temperature and 

FRICTIONAL HEAT GENERATION RATE 

VS 

LUBRICANT FLOW RATE 


CONSTANT FLOW RATE -226.79 KG/ HR 

(500 LBM/ H« ) 
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FIG.30 INNER RACE LUBRICANT TEMPERATURE, 
LUBRICANT MEAN TEMPERATURE AND 
CAGE SLIP PERCENT 

VS 

RADIAL LOAD 



e 

K 
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SEMi-LOGARITHMIC 359-71 
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Contract 


TABLE I 


ENCLOSURE 1 


NAS3-15700, Amendment No. 1 ■ 

; . ■ • • ^ .* *•••• 

ATTACIB-IENT I. ' ' 

(Revised 11-72) 

- BEARING SPECIFICATIONS AND OPERATING CONDITIONS 


Bore - 4.750 in. V . ■; 

Width - 1.875 in. . ‘ ' ^ , ’ 

o.D. - 8.125 in. •■■■■•. 

Roller length - (to be optimized) ' 

■. Roller large end diameter - (to be optimized) 

Roller end radius - (to be optimized) 

‘ Contact angle (cup) - (to be optimized) 

Mean Pitch diameter - (to be optimized) ' • ’ 

' Raceway surface finishes - 6 /^in« rms 
Roller surface finish - 4 /c; in. rms 

Large end flange s\irf ace finish - • 6 /< in . rms 
Number of rollers- - l4 (approximately) 

Inner ring, outer ring and roller steel - CVM SAE 4320 
• Shaft speed - 12,500 rpm • 

■ • Radial load - 2000-6000 lbs. ‘ • r. ! ' 

Thrust load - 2000-12,000 lbs. 

Lubricant temperature - 195^ 
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APPENDIX A 

ROLLER DIMENSIONS FOR A SPECIFIED 
BEARING SECTION AND VARIABLE CONTACT ANGLE 
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NOMENCLATURE 


Symbol 

A 

APEX 

B 

C 

D 

Dm 

Ds 

E 

em 

HT 

I.D. 

L 

O.D. 

T 

W 

Z 

otl 

O^o 


Description 

Minimum radial inner ring section height; ra (in,) 
Bearing APEX length; m (in.) 

Minimum radial outer ring section height; m (in.) 
Minimum inner ring flange width m (in.) 

Roller large end diameter; ra (in.) 

Roller mean diameter; ra (in.) 

Roller small end diameter 

Minimum distance between the roller large end/outer 
raceway contact point and the outer ring side face 

Bearing mean pitch diameter; ra (in.) 

Bearing radial section height; ra (in.) 

Inner ring inside diameter m (in.) 

Roller total length; ra (in.) 

Outer ring outside diameter m (in.) 

Bearing section limits; ra (in.) 

Bearing total width ra (in.) 

Total number of rollers 

Inner raceway contact angle depression radians 

Outer raceway contact angle degrees or radians 

Bearing mean angle degrees or radians 

One half of the roller included angle; degrees or 
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Computer program AE71Y001 requires full definition of bearing 
internal geometry as input. In order to use the program to 
evaluate bearing performance at a given contact angle it is 
required to establish a unique set of roller dimensions which 
will provide maximum load carrying capacity and which will fit 
within the bearing radial section. The following analysis was 
used to determine roller geometry for outer raceway contact angles 
of 12, 18, 20, 22, 24, 28°. 

The knowns in this analysis were the bearing boundary 
dimensions; O.D. = 8.125 in., I.D. = 4.750 in., total width, 
w = 1.875 in., and the outer raceway contact angle ( o<o ). In 
addition, referring to Fig. lA, dimensions A and B were both set 
equal to 0.2 inches, C was set equal to 15 percent of the total 
bearing width resulting in C = 0.282 inches and E was set to 
0.025 inches. The above dimensions are consistent with good 
design practice. 

By assuming values for half the roller included angle (“V> ), 
(see Fig. 2A) and the roller length to large end diameter ratio 
(L/D); the roller large and small end diameters plus the roller 
apex length can be found. 

^ ~ o — 

(lA) 

Since A equals B: 


D 


/ CO ^ 'Cx 

- -b (tAiJ t -TA.vJ 
"TAN V ^ ' 

D 

= X) - 2_LtamV 


(2A) 


(2A) 


C3A) 


Apex -- , 

2. S* fO V (4A ) 
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By defining different pairs of L/D and “V , a large but 
finite set of different roller geometries was generated for each 
contact angle. 

From the sets of roller dimensions generated, it remained to 
pick the particular set of dimensions which would fit within the 
bearing section and which would yield the maximum dynamic 
capacity for the given contact angle. 

-j- - O . "D . ~ T . - 

£ (5A) 

The following inequalities are used to insure that the 
roller geometry is such that it remains within the available 
bearing section. 

tests on radial height 

+ L sifJ 'C ^ Ut-A-B ("^i) 


T 

2 


tests on axial length considering inner ring flanges 


cof>-u^ ' 


T tests on axial length considering the outer ring width 
o 




L Cos t, S.W-C-E. 


T, 


If all three conditions are met, then the roller size is 
compatible with the bearing section. 


From the original set of roller dimensions generated, many 
have failed to pass one or more of the above inequalities and 
are eliminated from consideration. However, many sets of roller 
dimensions remain. From the remaining rollers, that roller 
which provided maximum dynamic load carrying capacity was selected 
as the optimum roller design for the given contact angle. _ A 
simplified equation was used to calculate the capacity of each 
remaining roller. 
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Having determined the roller size, the maximum number of 
rollers which can be put into the bearing is determined by 
considering the minimum distance allowed between rollers. In 
common practice, this distance is equal to 20 percent of the 
mean roller diameter. This allows for adequate cage web 
thickness. Therefore, the number of rollers- z -may be. obtained . as : 




DM is the roller mean diameter, e^^j is the bearing mean pitch 
diameter, Z must be rounded down to the next lowest integer. 


Fig. 3A is a plot of the roller dimensions versus outer 
raceway contact angle which resulted from this analysis. 


It is interesting to note that for the roller dimensions 
selected at each of the contact angles only two of the three 
dimensions T1 , T2 and T3 are nearly totally utilized. For 
contact angles 20° and lower dimensions T1 and T2 were nearly 
fully utilized. For contact angles 22° and above dimensions 
T1 and T3 were those fully utilized. 


119 

research laboratory SKP industries. INC. 



AL73P010 



FIG. I A BEARING SECTION 
LIMITING DIMENSIONS 

AND 

ACCEPTANCE LIMITS T 
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FIG.3A ROLLER LENGTH, ROLLER INCLUDED ANGLE 
and roller large end DIAMETER 
VS 

CONTACT ANGLE 
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